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(S) Color reproduction method and apparatus. 

(57) Method and apparatus for color printing ac- 
cording to a printer table having high color 
smoothness for out-of-gamut colors. A color 
printer gamut edge is first determined and color 
primary values for colors within the printer 
gamut are calculated and inserted into the 
printer table. For transition colors, namely 
those colors outside the printer gamut edge but 
still within the printer table, color primary 
values are calculated by selecting a color from 
the printer gamut edge at a point which lies at a 
constant angle from the transition color in 
question. The same constant angle is used for 
each and every one of the transition colors. A 
border table is provided for colors outside of 
the printer table ; the color primary values for 
the border table are selected in the same man- 
ner as for the transition colors. According to the 
invention, it is possible to provide smooth color 
transitions for out-of-gamut colors, and in par- 
ticular it is possible to provide monotonically 
increasing lightness in both the printer table 
and the border table and thereby avoid undesir- 
able bands of darker colors within lighter col- 
ors. 
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RELATED APPLICATIONS 

This application is related to a copending European patent application (agenf s ref. 

2271730). entitled "Method And Apparatus For High Fidelity Color Reproduction", the contents of which are 
incorporated by reference as if set forth here in full. 

BACKGROUND OF THE INVENTION 

This application includes an appendix of computer program listings. 

A portion of the disclosure of this patent document contains material which is subject to copyright protec- 
tion. The copyright owner has no objection to the facsimile reproduction by anyone of the document or the pa- 
tent disclosure, as it appears in the Patent Office patent file or records, but otherwise reserves all copyright 
rights whatsoever. 

Field Of The Invention 

The present invention pertains to a method and apparatus for building and using look-up tables which de- 
termine the colors that a color printer prints in response to requests to print specif ic colors. The specific colors 
requested to be printed may include colors that are not printable by the printer. For those colors in particular, 
the colors in the look-up tables vary smoothly so that there are not large, discontinuous jumps in printed color 
for only a small change in requested color. In addition, lightness for those colors changes monotonically, that 
is, for requests to print increasingly lighter shade of colors the actual colors printed are also increasingly lighter 
without any undesirable dips towards darker colors. 

Description Of The Related Art 

Recently, as the availability of color monitors and color printers has increased, it is more and more com- 
monplace for a computer user to view a full color image on a color monitor and then to request a full color printout 
of that image on a color printer. 

However, color printers and color monitors form color images differently. Specifically, a color monitor is a 
light emitting device; colors are formed on color monitors by adding light from three color primaries, generally 
red, green and blue. Printed images, on the other hand, simply reflect ambient light; colors are perceived by 
the way ambient light is affected by three subtractive primaries, generally cyan, magenta and yellow (and some- 
times black). 

These processes are fundamentally different. As a result, the range of colors displayable on the monitor 
is different from the range of colors printable by a printer Figure 1 is the CIE 1931 chromaticity diagram showing 
the range (or "gamut") of colors displayable by a monitor (area "A") and the range (or "gamut") of colors printable 
by a printer (area "B"). As seen in Figure 1, the range of colors displayable on a monitor is generally greater 
than the range of colors printable by a printer. This is because a monitor is a light emitting device and is able 
to display colors with greater saturation. There are, however, some low saturation areas such as at area 10 
where a printed image, which is light-sub tractive, has greater color range than a monitor. 

Because of the difference between the ranges of printable and displayable colors, it has not heretofore 
been possible to print color images which are perceived as faithful reproductions of displayed color images. 
Specifically, it is simply not possible to print a color in areas like out-of-gamut area 11 which are outside the 
range "B" of printable colors. Accordingly, even though those colors may be seen on color monitors, they cannot 
be printed on a color printer 

In U.S. Patent 4,941,038, out-of-gamut colors were adjusted to colors on the exterior of the printer gamut 
which had the shortest vector distance to the out-of-gamut colors and which preserved the chroma and hue 
of the out-of-gamut colors. But because each out-of-gamut color was adjusted independently, out-of-gamut 
colors are printed with poor color smoothness, where small changes in commanded color can result in large 
changes in printed color In particular circumstances, poor color smoothness manifests itself as non-monotonic 
changes in luminance whereby the lightness of out-of-gamut colors does not increase smoothly and monoton- 
ically from dark to light but rather dips occasionally from light to dark. This results in a situation where colors 
which should merge smoothly and monotonically from dark to light in fact show undesirable bands of darkness. 

Figure 2 shows an example of non-smooth color transitions in the form of non-monotonic changes in light- 
ness. Figure 2 shows a cross-section along an arbitrary hue plane in the a* and b* axes of CIELAB space and 
projected alone the L* lightness axis. In Figure 2, 24 designates the edge of the color printer gamut. Within 
edge 24, colors are printable and outside edge 24 the colors are out-of-gamut, unprintable colors. 26 designates 
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a border within which out-of-gamut colors are mapped into printable colors on the edge 24 of the printer gamut. 

Each of rows 27 represents a table which gives CMY values to be printed in response to a command to print 

the color in the corresponding color space. For purposes of Illustration, Figure 2 only shows the luminance 

value L* which results when the CMY values are printed. 
5 Ordinarily, it is expected for luminance to increase smoothly in directions parallel to th L* axis. Thus, for 

column 28 which Is parallel to the L* axis, the luminance L* value increases smoothly from L* = 50 to 56. 
For some areas, however, like column 29, the luminance starts to rise from L* = 50 to 52, but then dips 

back to L* = 51 before rising to the final level of L* = 52, Thus, the change in luminance in the resulting printed 

colors is non-monotonic showing undesirable dips or bands of darkness. 
10 Non-monotonic changes In lightness is only one example of non-smooth color transitions for out-of-gamut 

colors. 

SUMMARY OF THE INVENTION 

15 It is an object of the present invention to address the foregoing difficulties and to provide printer tables 

whose colors vary smoothly in the out-of-gamut regions, and in particular, whose colors in the out-of-gamut 
region exhibit monotonic increases in lightness. 

According to this aspect of the invention, transition colors are selected from colors on the edge of the printer 
gamut which are at a constant angle from the transition color without regard for change in lightness. Because 

20 each of the transition colors are at a constant angle from the edge of the printer gamut, the situation shown 
in Figure 2 cannot arise. Accordingly, out-of-gamut colors vary smoothly and exhibit monotonic increase in light- 
ness for increasingly lighter out-of-gamut colors. 

This brief summary has been provided so that the nature of the invention may be understood quickly. A 
more complete understanding of the invention can be obtained by reference to the following detailed description 

25 of the preferred embodiment thereof in connection with the drawings which together form a complete part of 
the specification. 

BRIEF DESCRIPTION OF THE DRAWINGS 

30 Figure 1 is a chromaticity diagram showing how the gamut of colors printable on a printer is related to the 

gamut of colors displayable on a monitor. 

Figure 2 illustrates printer tables in a cross-section of CIELAB color space. 
Figure 3 is a block diagram of a printing apparatus according to the invention. 

Figure 4 Is a flow diagram showing how a printer driver in the Figure 3 apparatus selects CMYK values 
35 for a color printer. 

Figure 5 is a flow diagram for describing how the printer table and the border table are constructed. 
Figure 6 illustrates a typical division of CIELAB space into a printer table. 
Figure 7 shows how unconnected regions are removed from the printer table. 

Figure 8 shows radially concave regions in the printer gamut and Figure 9 illustrates how to remove con- 
40 cavities by making the printer table radially convex. 

Figure 10 shows how CMY values are selected for each cell in the transition region of the printer tables. 
Figure 11 illustrates the arrangement of border tables. 

Figure 12 shows the relation between the printer table and the border table in CIELAB space for an arbi- 
trary luminance value L*. 

45 Figure 13 shows the printer gamut at an arbitrary luminance value L* in the a* and b* axes. 

Figure 14 shows how hue angles are warped for printer tables and Figure 15 shows how hue angles are 
warped for border tables. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

50 

Figure 3 is a block diagram of a printing apparatus according to the invention. 

As shown in Figure 3, the printing apparatus includes a host CPU 20, a color monitor 30 and a color printer 
40. Host CPU 20 includes a processor 21 such as an 80286 microprocessor, a random access memory ("RAM") 
22 which provides working storage area to processor 21 , a read only memory ("ROM") 24 which provides static 
55 storage for processor 21, monitor driver 25 and a printer driver 26. Host CPU 20 is accessed by an operator 
via keyboard 27 which is connected through an interface 29 to processor 21 . Using the keyboard, an operator 
can cause processor 21 to execute stored program instructions which cause color images to be displayed on 
monitor 30 and which cause corresponding color images to be printed on color printer 40. 
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Other peripheral devices, such as disk drives, tape drives, color video interfaces, color scanner interfaces, 
etc.. nnay be provided for host CPU 20 but those other devices are not shown in th interest of simplicity. In 
cooperation with the stored program instructions executed by processor 21 , such devices permit, for example, 
a color image to be scanned into RAM 22 and displayed on monitor 30, the colors in the image to be manipu- 
5 lated, and the resulting image to be printed on printer 40, 

In accordance with stored program instructions, processor 21 derives a color image for display on monitor 
30. Processor 21 provides the color image to monitor driver 25 which in turn derives RGB values for each pixel 
in monitor 30. The RGB values are provided via interface 31 to the monitor 30 where those values are displayed. 
Upon request, processor 21 also feeds a color image to printer driver 26 for printing by color printer 40. 

10 Printer driver 26 derives CMY values for each pixel of the color image based on the color values provided from 
processor 21 . The CMY values are determined in accordance with either a printer table 26a or a border table 
26b. The printer table 26a is a table which provides CMY values for all colors that are printable by printer 40. 
The border table 26b is a table which provides suitable CMY values for colors that are not printable by printer 
40. The printer table may also include CMY values for some unprintable colors so as to smooth the transition 

15 from printable to unprintable colors. In addition, a black (hereinafter "K") value may also be derived. The CMYK 
values are fed via interface 41 to printer 40 where they are stored in bit map memory 42 within printer 40. The 
bit map memory 42 may store a full bit map image of the printed image or it may store only a band or partial 
bit map image. When sufficient color data is stored in bit map memory 42. a color printer head 44 reciprocates 
across a platen adjacent to a sheet of paper In the present embodiment, print head 44 includes 32 ink jet noz- 

20 zles arranged in a four column by eight row pattern. The nozzles in the first column all eject droplets of cyan 
ink; the nozzles in the second column all eject droplets of magenta ink; the nozzles in the third column all eject 
droplets of yellow ink; and the nozzles in the fourth column all eject droplets of black ink. The nozzles are con- 
trolled independently in accordance with the color data in bit map memory 42 such that in one reciprocation 
of print head 44 across the platen, eight rows of pixels are printed. 

25 Figure 4 is a flow diagram showing how printer driver 26 selects CMYK values from the color data provided 

by processor 21. In step S401. printer driver 26 receives RGB values for a location (x,y) in bit map memory 
42. In step S402. printer driver 26 derives device independent color coordinates from the RGB value. Prefer- 
ably, the device independent coordinates are CIELAB coordinates. This is because the coordinates in CIELAB 
space are perceptually uniform such that equal-sized intervals anywhere in CIELAB space correspond to 

30 equal-sized changes in perceived color. Moreover, since CIELAB space can be viewed in cylindrical coordin- 
ates in terms of hue, saturation and luminance, it is an intuitive space which is amenable to defining gamut 
maps. 

In step S403, the luminance coordinate is compressed at the extremes of the L* axis in CIELAB space. 
Compression step S403 may be performed explicitly by mathematical manipulation of the L* value from step 

35 S402 or implicitly by storing modified CMY values in the printer table and the border table. If performed im- 
plicitly, which is preferable in some instances, then both the printer table and the border table store pre-com- 
pressed values. More particularly, the printer table and border table can be arranged so that the values stored 
at, for example, luminance L* = 99 actually correspond to a luminance of L* = 94. Likewise, values-stored at, 
say. luminance L* = 7 actually correspond to a luminance of L* = 26. Values in the center of the luminance 

40 range, from for example L* = 38 through 90, remain unmodified. This arrangement results in luminance com- 
pression without the need for explicit compression. 

While compression step S403 is optional, it is nevertheless preferable to perform since it ensures that col- 
ors at extreme values of luminance are printed with perceptible changes in luminance. More particularly, be- 
cause monitor 30 displays colors with light emitting elements, it can display colors with higher values of lumi- 

45 nance than those achievable by printer 40. whose highest value of luminance is limited by the whiteness of 
the paper upon which the color image is formed. Further still, since monitor 30 can completely turn off its light 
emitting elements, it can display colors with lower values of luminance than those printable by printer 40. since 
even black ink reflects some ambient light Accordingly, to ensure that some color is printed, even at the highest 
and lowest luminance values, it is preferable to compress the luminance values determined in step S402 into 

50 a range that is printable by printer 40, 

In step S404. the L*. a* and b" coordinates derived in steps S402 and S403 are inspected to determine 
whether they fell within the range covered by printer table 26a. If the L*. a*, b* coordinates are within the range 
covered by printer table 26a, th n flow advances to step S405 which looks up the corresponding CMY values 
in printer table 26a at location L*. a*, b* (actually, the nearest location since only discrete values of L*. a* and 

55 b* are stored). On the other hand, if the L*, a*, b" coordinates are not within printer table 26a. then flow ad- 
vances to step S406 in which the hue angle 9 is d rived from the a* and b* values according to the following 
formula: 

9 = arctan (b * /a • ) 
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Flow then advances to step S407 which looks up corresponding CMY values in border table 26b at the 
nearest location which corresponds to the lunnlnance L* and the hue angle derived in st p S406. 

in either event, flow then advances to step S408 in which the CMY valu s are stored in bit map memory 
42 at location (x,y). If desired, the CMY values may be modified before storage, for example, by interpolation, 
5 so as to accommodate the difference between the actual L*, a*, b* values stored in the tables and the desired 
values calculated above. 

In step S409, printer driver 26 determines whether the bit map memory has been completed. If the bit map 
memory has not been completed, then flow returns to step S401 in which the next RGB value Is received for 
the next location {x,y) in bit map memory. On the other hand, if the bit map memory has been completed, or 

10 if a sufficient area of the bit map memory has been completed (such as an eight row long band corresponding 
to the eight rows of ink jet nozzles in head 44), then flow advances to step S410 where gamma correction is 
performed. Gamma correction adjusts the CMY values in bit map memory so as to achieve a uniform distrib- 
ution of luminance. In step S411, undercolor removal is performed so as to derive the black value for location 
(x,y) in bit map memory. Undercolor removal in the present embodiment may be performed by the simple ex- 

15 pedient of selecting the minimum value of CMY and assigning that value to the black value. Then, each of the 
CMY values is adjusted by subtracting the black value from it. 

The order of steps S410 and S411 is not critical and those steps may be switched, for example to accom- 
modate a particular color printing technique such as continuous tone, dither techniques or error diffusion. 
In step S412, color printing is initiated using the resulting CMYK values. 

20 Figure 5 is a flow diagram for describing how the printer table 26a and the border table 26b are constructed. 

The flow procedures shown in Figure 5 need only be performed once for each printer, or onc^ whenever it is 
desired to re-calibrate each printer. More preferably, the flow prociedures of Figure 5 are performed only once 
for a class of printers, such as printers of the same model number, and are provided in software to an operator 
as part of a factory calibration of the printer. 

25 In step S501 , the gamut or range of colors printable by printer 40 is measured. Preferably, this is achieved 

by printing a very large subset, or a complete set, of all colors printable by printer 40. For example, in the printer 
used in the present embodiment, each of the CMY and K values may be printed in 65 gradations ranging nu- 
merically from 0 through 64, A subset of about one quarter of those values, for each color, are printed. Thus, 
for example, 17 C values are printed, namely numerical values 0, 4, 8, 12, ... 64, and 17 M values are printed, 

30 and 17 Y values are printed. All possible combinations of those 17 CMY values are printed, yielding 17x 17x 
17 = 4,91 3 color patches. 

In addition to the foregoing hued colors, all possible gray values, in this case 48 additional gray values 
over the 17 already printed, are also printed. 

With the foregoing sampling of the printer gamut, it will be seen that pure gray colors are printed together 
35 with hued colors. Whatever method of sampling is used, this property of pure gray printing should be preserved 
since proper gray reproduction is a desirable property in color reproduction. 

The color of each of the 4,913 color patches and 48 additional gray patches is measured in a device in- 
dependent color space such as the aforementioned CIELAB color space. Thus, at the end of step S501, for 
each of the 4, 91 3 + 48 = 4,961 unique CMY color combinations, L*, a* and b* coordinates are measured thereby 
40 defining the printer gamut. 

Step S502 derives mathematically smooth functions which map the CIELAB coordinates into CMY coor- 
dinates. In the present embodiment, a cubic least squares fit from CIELAB space into CMY space was chosen. 
That is, using well-known least squares fitting techniques, coefficients Cq through c^g, mo through mig, and yo 
through y^Q were derived to give the best fit, in the least squared sense, to the gamut measured in step S501: 

45 

C = Co + c^L* + Cja* + c-,b* + c^L*^ + Cja^^ + c^b*^ 
+ C7L*a* + CfX»*h* + C9a*b* + c,oL*^ + Ci,a*^ + 
^ c^b*' + c,3L*2a* + c,4L*a*2 + CjjL^^b* + c,<iL*b*^ 

+ c,7a*2b* + c,ga*b*^ + c,9L*a*b* (1) 

55 
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H- m7L*a* + ingL*b* + ii^a^b* + m,oL*^ + m,,a*^ + 
in,2b*^ + in,3L*^a* + in,4L*a*^ + m,5L*^b* + xn,^L*b*^ 
+ ro,7a*2b* + in,8a*b*^ + in,9L*a*b* (2) 



Y = Yo + YjL* + y^a* + yjb* + y^l.*^ + ysa*^ + y^b*^ 
+ yyL^a* + ygL*b* + y9a*b* + YioL*^ + Yii^*^ + 
y,2b*3 + y|3L*^a* + y|4L*a*2 + y,5L*^b* + y,^L*b*^ 
+ y,7a*^b* + Yiga^b*^ + y,gL*a*b* (3) 

In step S502, any mathematical function which fits the measurements taken in step S501 from the device 
independent coordinate space to CMY coordinate space may be used. Preferably, however, the mapping func- 
tion includes smoothing so as to eliminate measurement irregularities that may have been encountered in step 
S501. 

It may, in addition, be preferable to weight some of t he points measured in step S501 prior to deriving map- 
ping in step S502. For example, proper skin tone color reproduction is an important property of color printers. 
Accordingly, it may be desirable, in some circumstances, to weight colors in the area of skin tone colors more 
heavily than other colors. 

In step S503, the device independent space, namely CIELAB space, is divided into equally sized intervals, 
one of the intervals including the L* axis such as by being centered at the L* axis, thereby providing a blank 
printer table. The size of the printer table preferably includes both the printer gamut as well as the gamut of a 
typical color monitor. For example, referring to Figure 1, the printer table preferably includes the color area 
indicated generally at 12. The size of the intervals in the printer table should be made as small as possible 
giving due consideration to storage limitations for the printer table. Thus, for example, it has been found that 
fine luminance gradations are more important than fine hue and saturation gradations. It has been determined 
that dividing the luminance axis into intervals of AL* = 1 (luminance L* ranges from 0 through 100) provides 
adequate luminance gradation. On the other hand, such fine gradations are not ordinarily needed in hue, and 
therefore Aa* = Ab* = 3 has been found to yield adequate hue gradations (a* and b* range from about -100 
through +100 near the center (L* = 50) of the luminance axis). 

In addition to the foregoing considerations, it should also be observed that the printer gamut is not the 
same for each luminance value. Specifically, the gamut is relatively smaller at luminance extremes and rela- 
tively larger at the center of the luminance axis. Figure 6 illustrates a typical division of CIELAB space into a 
printer table, although all luminance and hue gradations have not been shown to simplify the presentation. At 
relatively low luminance values such as L* = 10, a relatively small rectangular grid in the a* and b* axes is ade- 
quate to map the printer gamut Similarly, at relatively high luminance values, such as L* = 90, a relatively small 
rectangular grid in the a* and b* axes is also adequate to store the printer gamut. However, at intermediate 
luminance values, such as that at L* = 50, a relatively larger rectangular grid in the a* and b* axes is required 
to map the printer gamut. 

As further shown in Figure 6, the rectangular grid at each luminance level includes the L" axis (in Figure 
6 it is centered on the L* axis). That is, there is a cell in the rectangular grid that corresponds exactly to a* = 
b* = 0. That central point, namely a* = b* = 0, corresponds to a pure gray color which, as mentioned above, is 
desirably reproduced as a pure gray color for proper color reproduction. 

In practice, it is also preferable to include more colors in the printer table than are in the printer gamut, 
and most preferably also to include the colors that are found in a typical monitor's gamut This permits the 
printer table to include transition values which smooth the transition from colors at the edge of the printer gamut 
to colors in the border table and which preserve color differ ntiation in areas outside the printer gamut 

In step S504, C, M and Y values are inserted into the printer table around the L* axis using the mapping 
functions derived in step S502. For digital color printing as opposed to continuous tone printing, fractional C, 
M and Y values are truncated or rounded to integral values. The entire rectangular grid at each luminance level 
is not completely filled, but rather only those points known to be within the printer gamut In addition, the points 
exactly on the L* axis, namely those points at a* = b* = 0, are also not mapped by the mapping functions. Rather, 
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CMY values for those paints are inserted in st p S505 by determining printer grays measured in step S501 
with corresponding L* values. This ensures that the smoothing introduced by the mapping functions does not 
introduce hue into pure gray values. 

In steps S506 and S507, the CMY values in the printer table are adjusted for unprintable colors. Unprintable 
5 colors may arise because of artifacts introduced by the mapping function selected in step S502. For example, 
the mapping function used may give rise to false regions in the printer table, such as region 45 in Figure 7, 
which is not within the printer gamut 46. These artifacts are removed in step S506 by removing all regions that 
are not connected to the region around the L* axis. 

Unprintable colors may also arise from the situation shown in Figure 8 in which 47 designates the edge 

10 of the printer gamut for an arbitrary luminance value L*. The printer gamut shown there is not radially convex 
because each and every radial line from the axis does not intersect edge 47 at one and only one point. In 
particular, radial line 48 intersects edge 47 at three points 49a, 49b and 49c. The region between points 49a 
and 49b is a radial concavity and can cause the generation of inappropriate CMY values in the printer table. 
Accordingly, in step S507, the values in the printer table are adjusted to make them radially convex. 

15 Figure 9 illustrates this process. Figure 9 shows the rectangular grid in the a* and b* axis for an arbitrary 

luminance value L*. Cells 51 through 55 all include printable values within the printer gamut. However, cell 59 
is a radially concave cell because the radial line at 9 crosses two celts in the printer gamut (cells 53 and 56). 
Accordingly, a CMY value is assigned to cell 59 to make the table radially convex. The value is selected by 
preserving as much as possible the hue of the color (angle 9 in Figure 9) and by selecting a saturation value 

20 which is nearest to that desired. Thus, in Figure 9, the value of cell 59 could be assigned a value near to the 
value of one of cells 51 through 55 in dependence upon which of cells 51 through 55 have the closer values 
of hue and saturation. In Figure 9, C = 1 , M = 1 8 and Y = 14 has been selected. 

Step S508 derives transition colors for each printer table, and Figure 10 shows how the CMY values for 
the transition colors are selected. In Figure 10, 60 is the edge of the printer gamut and 61 is the edge of the 

25 printer tables which, as mentioned above, corresponds roughly to the edge of a typical color monitor. For each 
transition color in the region between edges 60 and 61 , the color on the edge 60 of the printer gamut which 
lies at a constant anglefrom the transition color is selected for the transition color. Any change in luminance is 
allowed, that is, the change in luminance is not limited to an arbitrary threshold. For example, for transition 
color 62 point 64 on the edge 60 of the printer gamut is selected because point 64 lies at an angle a from point 

30 62. Likewise, for each point in the transition region, such as point 65, the corresponding point on the edge 60 
the printer gamut which lies at a constant angle a is selected, here point 66. For transition colors above the 
maximum saturated point 67 of the printer gamut, the angle a extends downwardly; conversely, for transition 
colors below the maximum saturated point 67, such as color 69, the angle a extends upwardly. For transition 
colors in the wedge 70 subtended by the angle 2a from the maximum saturated point 67, the maximum satu- 

35 rated point 67 is selected. Thus, alt colors in the wedge 70 are mapped to color 67. 

To ensure maximal color smoothness for the transition colors between edges 60 and 61 , the color satur- 
ation of the printer gamut edge 60 is preferably inspected before constant a angle extension to ensure that 
the color saturation increases monotonically from minimum saturation to the maximum saturation point and 
then decreases monotonically to minimum saturation. If non-monotonicity is found, then the color saturation 

40 at the printer gamut edge is adjusted before constant angle extension to remove the non-monotonicity. 

An angle a = 15** has been found to yield satisfactory results by allowing a pleasing increase in color sat- 
uration without unreasonably large changes in lightness. Other values of a, for example between 10"* and 20°, 
may also be used. 

Step S509 derives the border table 26b. Whereas printer table 26a was arranged as a rectangular grid in 
45 the a*, b* axes for each luminance value, the border table is arranged in wheels with one wheel for each of 
the luminance values of the printer table. Thus, as shown in Figure 11, a wheel-shaped border table is provided 
for each of the luminance values for which a printer table exists, which provides one border table in correspon- 
dence with each of the printer tables. The border tables contain plural cells which are accessed based on hue 
angle 9 calculated as a function of the a* and b* coordinates, as follows: 
50 9 = arctan(b*/a*) 

Figure 12 shows the correspondence of border tables and printer tables. Whereas printer table 26a is a 
rectangular grid for an arbitrary luminance value L*, border table 26b is arranged as a wheel centered at a* = 
b* = 0. Individual cells in the border table are accessed by the angle 9 in the a* and b* axes which also cor- 
responds to hue. Experimentally, it has been found that 1 degree increments, resulting in 360 cells in each 
55 border table, provides adequate gradation of hue, but this can be adjusted as described below in connection 
with Figure 13. 

Step S51 0 selects CMY values for each cell in the border tables in the same manner as for selecting CMY 
values for transition colors in the printer table. Thus, border table colors are selected from the printer edge color 
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which lies at a constant angle a (which is the same a as used for the printer table) from the border table. As 
befor , a extends upwardly or downwardly depending on whether the border table color is below or above the 
maximum printer saturation point, and is limited to th maximum printer saturation point when the border color 
table falls in wedge 70. 

In step S51 0, the border table values are inspected to ensure that the saturation of the colors in the border 
table changes smoothly. This situation is illustrated in Figure 13 which shows the printer gamut 70 at an arbi- 
trary luminance value L* in the a* and b* axes. As described above, colors such as color 71 which lie outside 
of the printer gamut are mapped to border color 72 of the printer gamut by preserving hue angle. Particularly 
in areas like area 74 where the printable saturation changes rapidly, small changes In hue can cause rapid 
changes in saturation in the border table. For example, as the hue angle changes from Oi to Oa. there are rel- 
atively large changes in saturation for only a small change in hue. Such a large change in saturation, when 
printed, looks unnatural. 

To avoid this unnatural look, the border table is increased in size until there are. sufficiently fine increments 
to ensure that saturation changes smoothly. If the border table size is increased, then the calculations in step 
S509 are repeated to fill in CMY values for the new border table. 

In step S511. the CMY values in the printer table are inspected and modified so as to ensure that they 
merge smoothly into the gray (L* axis). Specifically, at each discrete luminance level, the colors in the printer 
table that are dose to the L* axis are redetermined so as to ensure that they converge smoothly to gray. 

In step S512, the printer table is rectanguiarized. More specifically, until this step. CMY values have be n 
inserted into the printer table only in areas within the printer gamut 60 (steps S504 and S505) and in the tran- 
sition region 61 between the printer gamut and the border table (step S509). In step 8512, the remaining cells 
of the printer table such as cells like 69 in Figure 12 are filled out by calculating the hue angle for each blank 
cell that remains in the printer table and by inserting the border table color at that hue angle as illustrated at 
68 in Figure 12. 

In step S513, the hue angles in the printer table and in the border table are warped so as to compensate 
for the Abney effect. More specifically, the CMY values for out-of-gamut colors that are currently stored in the 
printer table and the border tables are all based on constant hue angle extension as described in connection 
with step 8508- However, for highly saturated colors, constant hue angle extension back to the printer gamut 
edge can result in a change in perceived hue, in accordance with the Abney effect. For example, constant hue 
angle extension changes the hue of a highly saturated (but unprintable) purplish-blue color to a less-saturated 
purple color on the printer gamut edge. 

To compensate for this effect, the hue angles in the printer table and in the border table are warped. More 
particularly, for both the printer table and the border table, the CMY values for one hue angle are transferred 
to another, different, hue angle so as to preserve the perceived hue of the printing color. This warping is illu- 
strated in Figure 14 and 15 for the blue/purple region of color space which lies at a hue angle 9 of between 9 
= 255^ to e = 333°. 

Figure 14 shows a printer table 80 before warping and the same printer table 81 after warping. The printer 
tables illustrated in Figure 14 are for arbitrary printer tables in the a*, b* plane and for an arbitrary luminance 
value L*, and it should be understood that the warping shown in Figure 14 is carried out for each of the printer 
table for the L* values selected in step S503. As shown in Figure 14, blue region 82a is stretched into warped 
blue region 82b. More particularly, each of the CMY values in region 82a is stretched into corresponding warp- 
ed regions in 82b. This ensures that the CMY values that are printed in response to a command to print a highly 
saturated out-of-gamut blue color yield a blue hued color rather than a purplish-biue color. For example, highly 
saturated out-of-gamut blue color 84, if printed according to the unwarped table 80 would yield a purplish-blue 
color while when printed according to warped table 81 yields a blue color. 

Additional hue angle warpings map the CMY values from region 85a into region 85b and map CMY values 
from region 86a into region 86b. The precise mappings are given as follows: 
For hue angles between 255° and 305° : 

warped_angle [255 + ang] = 255 + .5 * ang 
(where 0 < ang < 50°^ 

For hue angles between 305° and 309° : 

warped_angle [305 + ang] = 280 + 1.25 * ang 
(where 0 < ang < 

For hue angles between 309° and 333° : 

warped_angle [309 + ang] ~ 285 + 2 * ang 
(where 0 < ang < 24°) 

Thus, the unwarped region from 255° to 280° is warped by stretching into the region 255° to 305°, the un- 
warped region from 280° to 285° is warped by compressing into the region from 305° to 309°, and the unwarped 
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region from 285° to 333*' is warped by compressing Into the region from 309° to 333°, The warped regions re- 
main continuous, however, and they have the same end points as the unwarped regions (here 255° and 333°). 

Similar warpings are carried out in the red and cyan regions. For red, the warpings are: 
For hue angles between 10° and 40° : 

warped_angle [10 + ang] = 10 + 0.5 * ang 
(where 0 < ang < 3(P) 

For hue angles between 40° and 53° : 

warped_angle [40 + ang] = 25 + 1.25 * ang 
(where 0 < ang < 1T} 

For hue angles between 52° and 76°: 

warped_angle [52 + ang] = 40 + 1.5 * ang 
(where 0 < ang < 24°) 

Thus, the unwarped region from 1 0° to 25° is warped by stretching Into the region 1 0° to 40°, the unwarped 
region from 25° to 40° is warped by compressing into the region from 40° to 52°, and the unwarped region 
firom 40° to 76° is warped by compressing into the region from 52° to 76°. 

For cyan, the warpings are: 
For hue angles between 170° and 195° : 

warped_angle [170 + ang] = 170 + 2.0 * ang 
(where 0 < ang < 25°j 

For hue angles between 195° and 245° : 

warped_angle [195 + ang] = 220 + 0.5 * ang 
(where 0 < ang < 50°) 

Thus, the unwarped region from 170° to 220° is warped by compressing into the region from 170° to 195°, 
and the unwarped region from 220° to 245° is warped by stretching into the region from 195° to 245°. 

The same warping that is carried out on the printer tables is also carried out on the border tables as shown 
in Figure 1 5, Like before, this ensures that a request to print an out-of-gamut blue color, such as color 87, which 
would be printed as a purplish-blue hue if printed according to the unwarped border table, will result in a blue- 
hued color when printed according to the warped border table. 

With the warping described above, color smoothness in the printer table and in the border table is pre- 
served since both in-gamut and out-of-gamut colors are both warped. While it is possible to warp only out-of- 
gamut colors, this would result in an undesirable color discontinuity at the printer gamut edge. Moreover, even 
though the in-gamut colors are warped and hence are distorted, it has been confirmed experimentally that 
the amount of distortion is relatively insignificant because the hue angle rays are dose together for colors in 
the printer gamut and are further apart for higher saturation out-of-gamut colors. 

The above-described warping technique warps all color at the same hue angle equally regardless of the 
saturation of the colors. It is also possible to introduce a saturation-dependent warping factor, whereby more 
highly saturated colors are warped more than relatively less saturated colors. 

In step S514, colors in the yellow region of the printer and border tables are adjusted so as to widen the 
yellow region. More particularly, as shown in Figure 2, pure yellow colors for the printer fall into a very narrow 
range of the printer gamut that has been found to be difficult for users to find (the range of monitor yellow colors 
is greater). Because the range of pure yellow printer colors is so narrow, a user typically obtains a greenish- 
yellow rather than the desired pure yellow color. Thus, in step S514, the yellow region is widened. Conveniently, 
yellow widening is obtained through hue angle warping as follows: 
For hue angles between 87° and 91 ° : 

warped_angle [87 + ang] = 87 + 1.25 * ang 
(where 0 < ang < 4°) 

For hue angles between 91° and 97°: 

warped_angle [91 + ang] = 92 
(where 0 < ang < 6°) 

For hue angles between 97° and 112° : 

warped_angle [97 + ang] = 92 + 0.5 * ang 
(where 0 < ang < 15°) 

For hue angles between 112° and 132° : 

warped_angle [112 + ang] = 99.5 + 1.25 * ang 
(where 0 < ang < 20°) 

For hue angles between 132° and 147° : 

warped_angle [132 + ang] = 124.5 + 1.5 * ang 
(where 0< ang < 15°) 
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A computer program for automatically performing the foregoing steps S501 through S514 has been 
veloped and follows as an appendix to this description. 



MAKE -LOT {niake_lut:.c) 

Copyright (C) 1992 CANON INFORMATION SYSTEMS, INC. 
Copyright (C) 1992 CANON INC. 

This program uses the mathematical fit function from CIEIAB to CMY, 
which was derived from the print patches measurements. 
It builds the "LOT" file, which lists the CMYs for all I*AB intervals 
Together with the LOT it provides the '•lab_borders " file, which 
shows the rectangular border of the LOT gamut at each lightness leve 
It also builds the "border_colors " file, which lists the CMYs for th 
out of gamut colors at all angles and lightness levels. 
All these files are binary. 

-- Brigitte Ruetz 

***'^-Cc-et* ************************* **-tc*Kt*******^*Xt1eit**ir**'ttXe*** 



8 include "lut.h" 

extern convert (); 
extern cmy cmy_vector; 

static int startL, endL; 

static int inside^gamut [L^STEPS] (A_STEPS] [B_STEPSJ ; 
static int transition [L_STBPSJ [A_STEPS] CB_STEPS3 ; 
static int concave [L_STEPS] CA_STEPS3 [B__STEPS] ; 
static cmy LOT{L_STEPS] (A_STEPS] CB_STEPS] ; 
Static cmy border_color (I*_STEPS1 CANGliB_STEPS 3 ; 

static double interval_l tL_STEPSl , interval_a [A_STEPS3 , intezval__b [B_STEPS3 
static double saturation_l lANGLE_STEPS3 ; 
static double warped_angle CANGLE_STEPS3 ; 
static a_b rotated; 
static cmy averaged; 

int angle (a_index, b_index) /* computes the polar coordinates angle */ 

int a__index, b_index; 



{ 



int original_a, original_b, rounded_angle ; 
double fulltwo_angle, exact__angle ; 

original_a = a_index - CENTER_A; 
original_b ^ b_index - CENTER_B; 

fulltwo_angie = atan2pi ( (double) original_b, (double) original_a) 
if <fulltwo_angle >« 0) 

exact_angle « (ANGLE_STEPS/2) *fulltwo_ajigle; 
else 

exact_angle « (ANGLE_STEPS/2 ) * (f ulltwo_angle ♦ 2.0); 
rounded_angle « (int) (exact_angle + 0.5) ; 
if (rounded angle ANGLB_STEPS) 

rounded_angle O; 
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return rounded_angle; 

} 



double radius (ab_vector) /* computes the polar coordinates radius */ 

a b ab vector; 

{" " 

return sqrt (ab_vector.a*ab_vector.a + ab_v€ctor.b*ab_vector.b) ; 



rotate (angO, pointO) /* rotates a vector in the plane */ 

double aixgO ; 
a_b pointO; 

double CO, si; 

angO (2*ang0) /ANGLE^STEPS; 
CO « cospi(angO); " 
si « sinpi{angO); 

rotated, a = co*point0.a - si*pointO,b; 
rotated. b « si*point0.a + co*pointO,b; 

} 

average2(10. al, bl, a2, b2. round) /* computes the average of two colors 
int 10. al, bl, a2, b2, round; 
{ 

int angl, ang2; 

cmy vectorl, vector2; 

If (al 0 && al < A_STEPS && bl >« 0 && bl < B_STEPS 
35 && (inside_garaut [10] [al] [bl] == i | | transition [10] [al] [bl] 1)1 

vectorl.c - LOT(10] [al] [bl] ,c; 
vectorl.m = LUT[10] [al] [bl] .m; 
vectorl.y = LUTdO] [al] [bl] ,y; 

else 
40 I 

angX = angle (al, bl) ; 
vectorl.c = border color [10] [angl] .c; 
vectorl.m = border__color [10] [angl] .m; 
vectorl.y « border color [10] [angl] ,y; 

} 

45 if (a2 0 && a2 < A_STSPS b2 0 && b2 < B_STEPS 

(inside_gamut [10] (a21 [b2l i | | transition [10] [a2] [b2] «== 1)) 

vector2.c = LUT[10] [a2] [b2] ,c; 

vector2.m = LUTdOl [a2] [b2] .m; 

vector2,y = IjUT[10] [a2] [b2] ,y; 

50 ^ 



55 
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else 
{ 

ang2 = angle (a2, b2) ; 

vector2.c « boi:der_color flO] {ang2] .c; 
vector2-m « border_color (10] Cang2] ,m; 
^ vect:or2-y = border__color (101 Cang2] ,y; 

averaged, c = (unsigned char) ( (int) (0 .5* ( <int) vectorl. c 

+ (int) veccor2 .c) round*© . 5) ) ; 
averaged. m = (unsigned char) ((int) (0 .5* ( (int) vectorl, m 

+ (int) vector2 .m) + round*0.5)); 
averaged. y = (unsigned char) ( (int) (0 . 5* ( (int) vectorl.y 

+ (int) vector2,y) + round*0,5)); 

} 



average4(lk, a)c, bk, round) 

int Ik, ak, bk, round; 

{ 

if (ak CENTER__A - 2 && ak CENTER_A + 2 
&& bk >= CENTER_B - 2 bk <= CENTER_B + 2) 

averaged. c « LUTdk] [ak] [bk] .c; 
averaged. m « liUTUk] (akj [bk] .m; 
^ averaged. y = LUT[lk] [ak] (bkj .y; 

else 
{ 

xf (concave [Ik, ak - 1, bk] == 0 && concave [Ik, ak + 1, bk] 0) 
average2(lk, ak - i, bk, ak + i, bk, ro\md) ; 

else 

average2(lk, ak, bk - 1, ak, bk + 1, round) ; 



} 



void fit (Id, ad, bd, ci_ptr, mijptr, yi_ptr) /* confutes the CMY fits 

double Id, ad, bd; 

int *ci_ptr, *mi_jptr, *yi_ptr; 

int ck, mk, yk; 

double a2 , a3 , ab , a2b , b2 , ab2 , b3 , al , a2 1 , bl , abl , b2 1 , 12 , al2 , bl2 , 13 ; 
double cd, md, yd; 

a2 =t ad*ad; 
a3 = a2*ad; 
ab = ad*bd; 
a2b = a2*bd; 
b2 = bd*bd; 
ab2 = ad*b2; 
b3 = b2*bd; 
al = ad* Id; 
a21 = a2*ld; 
bl = bd*ld; 
abl « ad*bd*ld; 
b21 - b2*ld; 
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12 » ld*ld; 
al2 = ad*l2; 
bl2 = bd*12; 

13 = 12*ld; 

cd « FITC(ad,a2,a3.bd,ab,a2b,b2,ab2,b3,ld,al,a21,bl,aJDl,b2l,12,al2,bX2,l3) 
md = FITM(ad,a2,a3,bd,ab,a2b,b2,ab2,b3,ld,al,a21,bl,abl,b21,12,al2,bl2,13) 
yd « FITy(ad,a2,a3,bd,ab,a2b,b2,ab2,b3,ld,al,a21,bl,abl,b21,l2,al2,bl2,13) 

--^ '-^ >= 0) 

(int) (cd + 0.5) : 



30 



if (cd 

ck = 
else 

ck = 
if (md 

mk « 
else 

mk « 
if (yd 

yk = 
else 

yk - 
if (ck 

ck » 
if (mk 

mk = 
if {yk 

yk = 64 
if (ck == 

ck « 0; 
if (mk 

mk = 
if (yk 

yk - 0 
*ci_jptr 
*mi_ptr - mk; 
*yi_jptr = yk; 



(int) (cd 
>= 0) 
(int) (md 

(int) (md 
>= 0) 
(int) (yd 

(int) (yd 
== 65) 
64 



64 



65) 
65) 

-1) 



-1) 



0; 



-1) 
ck; 



- 0.5) 
+ 0.5) 

- 0.5) 
+ O.S) 



- O.S) 



35 



40 



/* intersects with printer gamut */ 



intersection ( 10 , angO , point 0 , 
updown, l_ptr, a_ptr, b_ptr) 
int angO, updown; 

double 10, *ljptr, *a_ptr, *b_ptr; 
a_b pointO; 

{ 

int Ik, ak, bk, ah, bh, al, bl, ck, mk, yk, count; 
int reached, final__ELh, final_bh, f inal_outside; 
double Is, 11, radO, radl, plane_ratio, warped_angO; 
double a^width, b^width, final_l, final_a, final_b, l_step; 
a_b place, placel; 



45 



50 



Is = saturation_l (angO) ; 
radO = radius (pointO) ; 

al = (int) ( (pointO -a - A^MIN) /DELTA_A) ; 
bl = (int) ( (pointO. b - B_MIN) /DEIiTA_B) ; 
l_step = niax2(0.0, 0.5* (10 - dark)); 
wa rped_angO = wa rped_angl e [ ang 0 ] ; 

if ( (warped_angO < ANGLE_STEPS/8 | | warped_angO >= 
(warped_angO >= (3*ANGLE_STEPS ) /8 && warped_angO 



( 7 *ANGI*E_STEPS ) /8 ) | | 
< ( 5 *ANGLE_STEPS ) / 8 ) ) 



55 



13 



BP 0 592 146 A2 



10 



{ ^ 

if^Sl^'^C^ER^Af ^ ' (2.0*warped_angO) /ANGLE_STEPS) ; 
{ 

for (ah = al; ah >= CENTER_A; ah--) 

place. a = interval_a [ahl ; 
place. b = plane_ratio*place.a; 
radl s= radius (place) ; 
if (radl > radO) 
continue; 

11 = 10 + min2(l step, -updown*SPACE_RATIO* CradO - radl) ) - 
rotate (angO - warped__angO , place); 
15 ^ (int) ( (rotated.a - A_MIN) /DELTA_A) ; 

bk = (int) ( (rotated, b - B__MIN) /DELTA_B) ; 
for (Ik = startL; Ik <« endL; lk++) 

+ 0,5* (interval 1 Uk + l] - interval 1 (IkJ ) > 1 
break; ~ _ * ^ ' 

} 

20 fitdl, place. a, place, b, &ck. Sink, &yk) ; 

if ((ck >= 0 && ck <= 64 && mk >= 0 && mk <« 64 
&& yk >i= 0 && yk <«= 64 
&& (inside__ganiut [Ik] [akj [bk] =« x 

I I (ak + 1 < A_STEPS 

&& inside^gamut [Ik] [ak + 1] [bk] == 1) 

II (ak + 1 < A_STEPS && bk + 1 < B_STEPS 

25 && inaide^gamut [Ik] [ak + X] [bk + l] x) 

I I (bk + X < B_STEPS 

&& inside_gainut [Ik] (ak] [bk + l] « x) 
I I (ak - X >- 0 && bk + X < B_STEPS 

&& inside_gaimit [Ik] (ak - XJ [bk + 1] x) 
I I (ak - X 0 

&& inside_gaimit [Ik] [ak - l] [bk] x) 
30 I I (ak - X >= 0 && bk - X 0 

&& inside^gamut [Ik] [ak - X] (bk - X] «= X) 
I i (bk - X >= 0 

&& inside_gaiirut Clk) [ak] [bk - x] X) 
I I (ak + X < A_STBPS && bk - X 0 

&& inside_gaimit[lk] [ak 4- i] [bk - i] =« x) ) ) 
II inside_ganiut [Ik] (ak) [bk] == i) 

35 { 

if {ck < 0 I I ck > 64 1 I mk < 0 1 | rak > 64 . 

I I ylc < 0 ! I yk > 64) 

f inal_outside = l; 
else 

f inal_outside = 0; 
final_l = 11; 
^ final_a = place. a; 

final_b - place. b; 
fiaal_cdi = ah; 
break; 

) 

else if (updown* (Is - 11) > o) 

^ f inal_outside =0; 

final_l =11; 
final_a = place, a; 
final_b = place.b; 
f inal_ah = ah ; 

50 
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break ; 
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} 

] 

else 
{ 



for (ah ^ al; ah <- CENTER A; ah++) 
{ 

place. a = interval_a [ah] ; 
place. b «= plane_ratio*place.a; 
radl = radius (place) ; 
if {radl > radO) 
continue; 

11 = 10 + niin2 (l_step, -updown*SPACE_RATIO* (radO 
rotate (angO - warped_angO, place); 
a3c = (int) ( (rotated. a - A_MIN) /DELTA_A) ; 
bk - (int) ( (rotated. b - B_MIN) /DELTA_B) ; 
for (Ik = startL; Ik <:= endli; lk++) 



radl) ) ; 



{ 



if (interval^l tlk] 
breaJc ; 



+ 0.5*(interval_l(lk + 1] - interval l[lk]) > 



} 

fitdl, place. a, place. b, &ck, &mk, &yk) ; 
if ( (ck 0 &5e ck 64 && oak >= 0 && mk 64 
yk >« 0 && yk <« 64 
&& (inside_gantut Clk] [ak] [bk] == l 
(ak + 1 < A_STEPS 
&& inside_gainut [Ik] [ak + 1] [bk] 1) 

(ak + 1 < A_STEPS && bk + 1 < B_STBPS 
&& inside_gainut [Ik] [ak + 1] (bk + 1] 1) 

(bk + 1 < B_STEPS 
&& inside_gamut [Ik] [ak] [bk + 1] 1) 

(ak - 1 0 && bk + 1 < B_STEPS 
&& inside_gamut [Ik] [ak - IJ [bk + 1] == 1) 

(ak - 1 >= 0 
&& insidejamut [Ik] [ak - 1] [bk] 1) 

(ak - 1 >« 0 && bk - 1 >= 0 
&& inside_ganiut [Ik] [ak - 1] [bk - 1] 1) 

(bk - 1 >« 0 
&& inside^gamut [Ik] [ak] [bk - l] l) 

(ak + 1 < A_STKPS && bk - 1 >= 0 
&& insidegamut [Ik] [ak + 1] [bk - 1] 1))) 
I 1 inside__gamut [IkJ [ak] [bk] i) 

if (ck < 0 I I ck > 64 I j mk < 0 | ] mk > 64 
I 1 yk < 0 I I yk > 64) 
f inal_outside = i; 
else 

f inaloutside = 0; 
final_l = 11; 
final_a = place, a; 
final_b = place. b ; 
final^ah = ah; 
break; 



} 

else If (updown*(ls - li) 

f inal_outside « 0; 
final_l ^11; 
£inal_a = place .a; 



> 0) 
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final_b = place. b; 
final_ah = ah; 
break ; 

} 

) 

f (f inal_outside -= i) 
reached ~ 0; 

if {warped_angO < ANGLE_STEPS/4 || warped_angO (3*ANGIiE_STEPS) /4) 

for <ah = final_ah; ah >= CEKTER_A; ah--) 

placel.a = int:erval_a [ah] ; 
placel.b « plauie_rat:io*placel.a; 
radl = radius (placel) ; 

11 10 + min2 (1 step, -updown*SPACE_RATIO* (radO - radl)}; 

fitdl, placel, a, placel.b, &c)c, &mk, &ylc) ; 

if ( (ck >« O && ck <« 64 && mk >= O && mk <= 64 

^ &£ yk >= 0 && yk <= 64) || (updown*(ls - ll) > 0)) 

reached 1; 

a^width = 0.5*(final_a - placel.a); 
place. a = placel.a + a_width; 
break; 

else 
{ 

for (ah = flnal_ah; ah CENTER_A; ah++) 

placel.a = interval a [ah] ; 
placel.b = plane_ratio* placel. a; 
radl ^ radius (placel) ; 

11 10 + miii2 (l_step, -updown*SPACE__RATXO* (radO - radl)); 

fitdl, placel.a, placel.b, &ck, &mk, &yk) ; 

if ((ck >= 0 && ck 64 && ink 0 && mk <= 64 

&& yk >« 0 && yk <« 64) | | (updo wn* (Is - ll) > 0)) 

reached = i; 

a_width = 0 .5* (f inal_a - placel.a); 
place, a = placel.a + a_widt:h; 
break ; 

if (reached =- 0) 

{ 

place. a = 0; 
place. b = 0; 

if (warped_angO < ANGLE_STEPS/4 || warped_angO >= {3*ANGLE_STEPS) /4 

a_width « max2 (f inal_a, incerval_a [CENTBR_A +11); 
else 

a_width = min2 (f inal_a, int:erval_a (CENTER_A - l] ) ; 

count = 0 ; 

while (count < 11) 

{ 

a width « 0.5*a width ; 
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place, b = plaiie_rat.io*place.a; 
radl =s radius (place) ; 

11 « 10 + miii2 (Instep, -updown*SPACE_RATIO* (radO - radl) ) ; 

if (place. a == 0 && place. b =^0) 

( 

if (10 < 30) 
{ 

final_l =11; 
final_a = 0; 
final b = 0; 

} 

place. a = place. a + a_width; 

else 
{ 

fit (11, place. a, place. b, &ck. &ink, fcyJc) ; 
if (ck >= 0 && ck <= 64 && mk 0 && mk <= S4 
yk >= 0 && yk <= 64) 

{ 

final_l = 11; 

final_a = place. a; 

final_b = place. b; 

pi ace. a - place. a + a_width; 

else if (updowii*(ls - ll) > 0) 
place. a « place. a + a_width; 
else 

place. a « place. a - a_width; 
count++; 

} 

} 

else 

if (warped_angO < ANGLERS TE PS/ 4 || warped^angO >« (3*ANGLE_STEPS) /4) 

for (aih = final ah; ah < A STEPS; ah++) 
{ 

placel.a ^ iatexval_a [ah] ; 
placel.b « plane_ratio*placel.a; 
radl = radius (placel) ; 
if (radl > radO) 

{ 

a_width = 0 . 5* (placel .a - f inal_a) ; 
place, a = placel.a - a^width; 
break; 

) 

11 = 10 + min2 (l_scep, -updown*SPACE_RATIO* (radO - radl)); 

fitdl, placel.a, placel.b, &ck, &£nk, &yk) ; 

if ( (ck < 0 II ck > 64 I I mk < 0 | | mk > 64 

^ 1 I yk < 0 I I yk > 64) |I (updown* (Is - 11) :> 0)) 

a_width ~ 0 .5* (placel.a - f inal_a) ; 
place. a = placel.a - a_viidth; 
break; 

) 

else 

{ 
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for (ah - final_ah; ah 0; ah--) 

placel.a = interval_a (ah] ; 
placel.b « plane__rat:io*placel-a; 
radl « radius (placel) ; 
if (radl > radO) 

a_width = 0,5* {placel.a - f inal_a) ; 
place. a = placel.a - a_widch; 
break; 

} 

11 = 10 -f min2 (l_scep, -updown*SPACE_RATIO* (r-adO - radl) ) ; 

fitdx, placel.a, placel.b. &clc, &nik, feyk) ; 

if < (ck < 0 II ck > 64 I 1 mk < 0 | | ink > 64 

^ I I yk < 0 I I yk > 64) || (updovin*(ls - 11) > 0)^ 

a__width =x 0.5* (placel.a - final_a) ; 
place. a « placel.a - a_width; 
break; 

} 

count = 0 ; 

while (count < ii) 

{ 

a_width » 0.5*a_widUi; 

place, b « plane_ratio*place.a; 

radl = radius (place) ; 

11 « 10 + min2(l_step, -updown*SPACE_RATXO* (radO - radl) ) ; 
fit(ll, place. a, place. b, &ck, 5ank, &yk) ; 
if (ck 0 && ck <« 64 && mk 0 && mk <= 64 
&& yk 0 && yk <» 64) 

final_l « 11; 

final_a = place .a; 

final_b = place. b; 

place. a = place. a + a_width; 

else if (updown* (Is - 11) > 0) 
place -a = place. a + a_width; 
else 

place. a = place. a - a_width; 
coxmt++; 

else 
{ 

if (warped_angO == ANGLE_STEPS/4 || warped_angO (3*ANGLE_STEPS) /4 ) 

plane_ratio = 0; 
else 

{ 

plane_ratio = 1 . 0/tanpi ( (2 . 0*warped_cUigO) /ANGLB_STEPS) ; 

if (bl > CENTER B) 
{ 

for (bh = bl; bh >= CENTER B; bh--) 

{ " 
place. b « interval_b [bh] ; 
place. a = plane_ratio*place . b ; 
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50 



radl ~ radius (place) ; 
if (radl > radO) 
continue; 

11 =^ 10 + min2(l_step, -updown*SPACE_RATIO* (radO 

rotate (angO - warped_angO, place) ; 

ak « (int) ( (rotated. a - A_MIN) /DELTA_A) ; 

bk = (int) ( (rotated. b - B_MIN) /DELTA^B) ; 

for (Ik = startL; Ik <= endLi; l)c++) 



radl) ) 



if (interval_l[lk] + 0. 5* (interval_l [Ik + IJ - interval_l [Ik] ) > l 
break; 



} 

fit (11, place. a, place. b, &ck, &mk, &yk) ; 

if ( (ck >= 0 && ck <« 64 && ink >= 0 && mk <= 64 

&& yk 0 £c& yk <= 64 

&& (inside^gamit [Ik] [ak] (bk] l 

I I (ak + 1 < A_STfiPS 

&& inside_gamut [Ik] (ak + l] [bk] 1) 

II (ak + 1 < A_STEPS bk + 1 < B_STEPS 
&& inside_gamut [Ik] [ak + i] (bk + i] i) 

II (bk + 1 < B_STEPS 

&& inside_gamut [Ik] [ak] [bk + 1] «« 1) 
11 (ak - 1 >= 0 && bk + 1 < B_STEPS 

&& inside_gamut [Ik] [cik - 1] (bk + 1] == l) 
I I (ak - 1 >« 0 

&& inside_gamut [Ik] [ak - 1] (b]c] 1) 

I I (ak - 1 >« 0 && bk - 1 >« 0 

&& inside_gainut[lk] (ak - 1] (bk - 1] == i) 

II (bk - 1 0 

&& inside_gainut [Ik] (ak] (bk - 1] 1) 
I I (ak + 1 < A_STEPS && bk - 1 0 

&& inside_gamut (Ik] (ak + l] (bk - 1] i) ) ) 
I I inside gamut [Ik] [ak] [bk] i) 

if (ck < 0 I I ck > 64 II mk < 0 | | znk > 64 

I I yk < 0 I I yk > 64) 

f inal_outside = 1; 
else 

f inaloutside = O; 
final_l « 11; 
final_a - place. a; 
final_b « place. b; 
final_bh = bh; 
break; 

} 

else if {updown* (Is - 11) > o) 



{ 



} 



} 



f inal__outside - O; 
final^l = 11; 
finala = place, a; 
final_b = place. b; 
final__bh = bh; 
break; 



) 

else 
{ 

for (bh = bl; bh <= 
{ 



CENTER B; bh++) 
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5 place. b = interval_b [bh] ; 

pi ace. a = plajie_rat io ♦place .b; 
radx = radius (place) ; 
if (radl > radO) 
continue; 

11 « 10 + niin2 (l_step, -updown*SPACE_RATIO* (radO - radl) ) ; 

rotate (angO - warped_angO, place); 

ak = (int) ( (rotated.a - A_MIN) /DEIjTA_A) ; 

bk « (int) { (rotated.b - B_MIN) /DEIjTA_B) ; 

for (Ik = startL; Ik <« endL; lk++) 

if (interval_i [Ik] + 0 . 5* (interval_l tlk + 1] - incerval_l (Ik] ) > 
break; 

15 ] 

fit (11, place, a, place. b, &ck, &mk, &yk) ; 
if ( (ck 0 && ck <= 64 && mk >= 0 && mk <= 64 
yk >= 0 &Sc yk <= 64 



20 



25 



30 



{ 



&& Cinside_gainuc [Ik] [ak] Cbk] == 1 
(ak + 1 < A__STEPS 
Sc& inside_gamut (Ik) £ak + 1] [bk] ~ 1) 

(ak + 1 < A_STEPS && bk + 1 < B_STEPS 
&& inside_gamuc [Ik] [ak + i] [bk + 1} x) 

(bk + 1 < B_STEPS 
&& inside^gamut [Ik] [ak] [bk + l] 1) 

I I (ak - 1 0 && bk + 1 < B_STEPS 
&& inside^gamut [Ik] [ak - 1] (bk + 1] 1) 

(ak - 1 >« 0 
&& inside_gamut [Ik] [ak - 1] [bk] 1) 

(ak - 1 >« 0 && bk - 1 >= O 
&& inside_gainut [Ik] [ak - 1) [bk - 1] 1) 

II (bk - 1 >« 0 
&& inside_gainut [Ik] [ak] [bk - X] 1) 

(ak + 1 < A STEPS && bk - 1 >= 0 

inside^gamut [Ik] [ak + 1] [bk - 1] =^ 1) ) ) 
inside_gainut [Ik] [ak] (bkj == X) 



if (ck < 0 i I ck > 64 j I mk < 0 i | mk > 64 
II yk < 0 I I yk > 64) 
f inal_outside = 1; 
35 else 

f inal_outside - 0; 
final_l « 11; 
final_a = place, a; 
final_b = place. b; 
final_bh » bli; 
break; 

40 y 

else if (updovm* (Is - 11) > 0) 

f inal_outsid€ = 0; 
final_l « 11; 
final a = place. a; 
45 final_b = place -b; 

final_bh - bti; 
break ; 

1 



50 



if (final outside == 1) 
{ 
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reached = 0; 

if (warped_angO < ANGIjE__STEPS/2) 

for (bh final_bh; bh >= CENTER_B; bh--) 

placel.b = int:erval_b [bh) ; 
placei.a = plane_rat:io*placei .b; 
radl = radius (placel) ; 

11 = 10 + min2 (l_step, -updown*SPACE_RATIO* (radO - radl)) 
fitdl, placei.a, placel.b, &ck, &mlc, &yk) ; 
if ( (clc >= 0 ck <= 64 && mk 0 mk <= 64 
^ ScSc yk >= 0 && yk <= 64) || (updown*(ls - 11) > 0)) 

reached = l ; 

b_width = 0.5*(final_b - placel.b); 
place. b = placel.b + b_width; 
break; 

else 

{ 

for (bh = final bh; bh <= CENTER B; bh++) 

{ " " 

placel.b = interval_b [bh] ; 
placei.a « plane_rat:io*placel.b; 
radl = radius (placel) ; 

11 = 10 + min2 (l_step, -updown*SPACE_RATIO* (radO - radl)) 
fit(ll, placei.a, placel.b, &ck, &mk, &yk) ; 
if ( (ck >= 0 && ck <= 64 && mk >= 0 && mk <= 64 
yk >= 0 && yk 64) | | (updown* (Is - 11) > O) 

reached = 1; 

b_width = 0.5*(final__b - placel.b); 
place, b = placel.b + b_widt:h; 
break; 

if (reached 0) 

{ 

place. a = 0; 
place. b = 0; 

if (warped^angO < ANGLE_STEPS/2) 

b_width = max2 (f inal__b, interval ^ [CENTER_B + 1]); 
else 

40 b width = inin2 (final b, interval b [CENTER B - 13); 

) ' ~ " 

count = 0 ; 

while (count < 11) 
{ 

b_width = 0.5*b_width; 
45 place. a - plane_ratio*place .b; 

radl = radius (place) ; 

11 = 10 + inin2 (l_step, - updown* S PACE_RATIO* (radO - radl)); 
if (place. a == 0 && place.b == 0) 

if (10 < 30) 
{ 

final 1 = 11; 



30 
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final_a « O; 
final b = 0; 
5 } " 

^ place. b = place. b + b__width; 

else 
{ 

fitdl, place-a, place.b, &ck, &mk, Siyk) ; 
10 if (ck >« 0 && ck <= 64 mk 0 && mk <= 64 

&& yk >= 0 && yk <= 64) 

final_l « 11; 

final_a « place .a; 

final_b = place. b; 

15 ^ place. b = place, b + b_width; 

else if (updown*(ls - 11) > 0) 
place. b « place -b + b^width; 
else 

place. b = place. b - b__width; 

20 } 

count ++; 

} 

} 

else 
{ 

25 if (warped_angO < ANGLE_STEPS/2) 

for (bh = final_bh; bh < B_STEPS; bh++) 

placel-b = interval_b (l>h] ; 
placei.a « plane_ratio*placel,b; 
3Q radl = radius (placel) ; 

if (radl > radO) 
{ 

b_widt:h - 0 .5* (placel .b - final_b) ; 
place. b ^ place l,b - b^width; 
break ; 

35 } 

11 = 10 + min2 {l_st:ep, -updown*SPACE_RATIO* (radO - radl) 
fitdl, placel. a, placel. b, &ck, &mk, &yk) ; 
if ( (ck < 0 II ck > 64 ! I mk < O | | mk > 64 

I I yk < 0 I I yk > 64) || (updown* (Is - 11) > 0)) 

40 b_width = 0 .5* (placel. b - final_b) ; 

place. b = placel, b - b_widch; 
break ; 

4S else 

for (bh = final bh; bh >= 0; bh- - ) 
{ 

placel. b = interval_b [bh] ; 
placel. a « plane_ratio*placel .b ; 
50 radl = radius (placel) ; 

if (radl > radO) 

{ 

b_width = 0,5* (placel. b - f inal_b) ; 
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place. b placei.b - b width; 
break; — ' 

} 



fit (11 Dla~? ^- -updown*SPACE_RATIOMradO - radl) ) ; 

I ^ ck > 64 I mk < 0 II mk > 64 

II < 0 II yk > 64) ( (updovrAids - ll) > 0)) 



b width = 0.5*(placel.b - final b) ; 
place. b = placei.b - b width- ~ 
^ break; ~ 

count - 0; 

while (count < ii) 

b_width = 0,5*b_width; 
place. a = plane_ratio*place,b; 
radl = radius (place) 



™^»2(l_step, -updown*SPACE RATIO* (radO - radl) ) • 
fxt 11, place.a, place. b, &ck, &ink, &yk) ; 
If (ck >=. 0 && ck <= 64 t& mk 0 && mk <« 64 



^ " ' - - u && mk <e 64 

^ yk >= 0 && yk <= 64) 

final_l = 11; 
final_a = place. a; 
25 final_b « place, b; 

^ place. b « place. b + b__width; 

else if {updown*{ls - li) > o) 
place. b - place. b + b width; 
else ~ 

30 place. b = place. b - b width; 

count++; ~ 

*1 _ptr = final_l; 
35 *a_ptr = final^a 

♦b_ptr = final_b 



40 



45 



main () 
{ 

PILE *lab borders_ptr, *lut_ptr. *border colors ptr; 
FILE *m_minBjptr, *m_maxB _ptr; ~ 
int min, max, first, last, index, sign; 
int ll, ai, bi, Ij, aj, b j , 11, aa, bb; 
mt ang, angl, ang2, ang_diff; 

S^^'u'i'^'wr' ^"-^ "^3' cc, mm, yy, kk; 

int h, h], hh; 

int final_ci, final_mi, final__yi, final ai, final bi- 
int lowL, mediumL, highL; " J.iicij._ui, 

int cyan, magenta, yellow; 
50 int start_ai[ANGLE_STEPS] , end_ai {ANGLE STEPS]- 

int start_bi [ANGLE_STEPS] , end_bi [ANGLE'stepS] • 
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int ^^o^e_a(8]| close^brs), far a [163 , far b(16]; 

int startML^TEPSl , endA[L_STBPSj , scartB[L STEPS], endB [ I. STEPS]; 

int minB[L_STEPS] (A_STEPS1 , maxB (L_STEPS] [A STEPS] ; ~ 

xnt new^minB CL_STEPS] [A_STEPS3 , new maxB[L STEPSl-fA STEPS] ; 

inc m minB[L_STEPSl [A__STEPS] , m ma^[L STEPS] [A STEPS] ; 

double left, right, down, up; ~ ~ ~ 

double saturation, new_saturatiou, factor- 

double 1, a, b, Im, In, an, bn, rad, radra' ratio, delta 1; 
double warped_ang, double_ang; ~ 
any hue [16], saturation_border_color (ANGLE_STEPS1 ; 

point, saturation_border (ANGbE_STEPS] ; 
a_b outer border [L_STEPS] [ANGLE STEPS] ; 

if ( (m_minB_ptr « f open ( "m_minB" , "r-)) NULL) /* opens files */ 



printf ( "cannot open m_minB\n"); 
exit (1) ; 
) 

20 ( <in_maxB_jptr = f open { »in_maxB" , «r")) NULIi) 

printf ( "cannot open m inaxB\n»*) ; 
^ exit(l); ~ 

if ( (lab_borders__ptr = f open ( lab_borders , "wb")) NUIiIj) 

25 printf ("cannot create lab borders \n «) ; 

exitCl); 

} 

if ({lut_ptr = fopendut, "wb")} ^« NULL) 

pr int f { "cannot create lut\n"); 
exit(l); 

30 } 

if ( (border_colors_ptr « f open {border_colors , "wb")) == NULL) 

print f ("cannot create border colors\n") ■ 
exit(l); 

) 

for (li = 0; li < L_STEPS? li++) /* ^eads files */ 

^ for (ai - 0; ax < A_STBPS; ai++) 

m_minBCli] [ai] = B_STEPS; 
ni_inaxB(li] [ai] = -1; 

for (li = 0? li < L_STEPS; li++) 

^ for (ai = m_startA[li] ; ai <« m_endA[li] ; ai++) 

fread{&m_minB[li] [ai] , sizeof (int) , 1, m rninB_ptr) ; 
^ fread(&m__inaxB[li] [ai] , sizeof (int), l, mlroaxBjptr) ; 

min = (int) ({dark - L_MIN) /DELTA_L) ; /* computes tile LABs of the 

max « (int)( (light - L__MXN) /DELTA_L) ; /* interval centers */ 

lowL min + 1 - extendL; 
mediuxnL « min + 1 + extendL; 

50 
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highL = max - 1 - (L_STEPS - max) ; 
for (li « lowL; li < mediumL; 1x4-+) 

interval_l Clil = Ij_MIN + (rain + 1 + (li - lowL)*0,5 + 0 , 2S) *DELTA_L; 
for (li mediumL; li < highL; li+-»-) 

intierval_l [li] = L_MIN + (li + 0 . 5) *DELTA_L; 
for (li =. highL; li <= L_STEPS; li++) 

iiiteirval_l [lil ^ L_MIN + (highL + (li - highL) *0. 5 + 0 - 25 ) *DELTA_L; 
for (ai = 0; ai < A_STEPS; ai++) 

interval_a[ai] = A_MIN + (ai + 0 . 5) ♦DELTA_A; 
for (bi = 0; bi < B_STEPS; bi++) 

interval__b{bi] ^ B_MIN + (bi + 0 . 5) *DELTA__B; 



for (ang = 0; ang < ANGLE_STEPS; ang++) /* computes warped angles 

warped_angle [ang] = ang; 
for (ang - 0; ang < 3; ang++) 

warped_angleCANGLS_STBPS/4 - ang) = yellow_ang - 1.25* (ang + l) ; 
ang_diff « yellow_ang - ANGLE_STEPS/4 ; 
for (ang =« 1; ang <= ang_diff + yellow_shi£t; ang++) 

warped_angle (ANGLE STEPS/4 + ang] - yellow_eUig; 
for (ang = X; ang <= ANGLE_STEPS/24 ; ang++) 

warped_angle [ANGLE_STEPS/4 + ang_diff + yellow shift + ang] 

« yellow__ajig + 0.5*ang; 
for (ang « 1; ang ANGLE STEPS/30; ang++) 

warped_angle (ANGLE_STEPS74 + ANGLE_STEPS/24 

+ ang_diff + yellowshift + ang] 

= yellow_ang + 0 . 5*ANGLE_STEPS/24 + 1.2S*ang; 
for (ang - 1; ang <= 2*(ang_diff + yellow shift) + 5; ang++) 

warped_angle [ANGLE_STEPS/4 + ANGLE_STEPS/24 + ANGLE_STEPS/30 

+ ang_diff + yellow_shift + ang] 

yellow_ang + 0 , 5*ANGLE_STEPS/24 + AKGLE_STEPS/24 

+ l.5*ang; 

angl = ANGLE_STEPS/36 ; 

for (ang « 0; ang < ANGLE_STEPS/l2 ; ang++) 

warped_angle [angl + ang] = angl + 0 .5*ang; 
angl = ANGLE_STEPS/24 + ANGLE_STEPS/36 ; 
ang2 « ANGLE__STEPS/12 + ANGLE_STEPS/36 ; 
for (ang 0; aiig < ANGLE_STEPS/18; ang++) 

warped_angle (ang2 + ang] ^ angl + 1.25*ang; 
angl = ANGLE_STEPS/8 + ANGLE_STEPS/72 ; 
ang2 ANGLE_STBPS/6 ; 

for (ang = 0; ang < ANGLE_STEPS/36 ; ang++) 
warped_cLngle (ang2 + ang] = angl + l.5*ang; 

angl = {5*ANGLE_STEPS) /8 + ANGLE_STEPS/18 ; 

for (ang = 1; ang <= ANGLE_STEPS/12 + ANGLE_STEPS/18 ; ang++) 

warped_angle [angl - ang] = angl - 0.5*ang; 
angl - ANGLE_STEPS/2 + ANGLE_STEPS/12 + ANGLE_STEPS/36 ; 
ang2 ^ ANGLE_STEPS/2 + ANGLE_STEPS/24 ; 

for (ang = 1; ang <= ANGLE_STEPS/ia + ANGLE_STEPS/72 ; ang + + ) 
warped_angle [ang2 - ang] = angl - 2,0*ang; . - 

angl = (3*ANGLE_STEPS) /4 - ANGLB_STEPS/36 ; 

for (ang = 0; ang < ANGLE_STEPS/12 + ANGLE_STEPS/36 ; ang++) 

warped_angle (atngl + ang] = angl + 0.5*ang; 
angl = (3*ANGLE_STEPS) /4 + ANGLE_STEPS/36 ; 
ang2 - (3*ANGLE_STEPS) /4 + ANGLE_STEPS/12 ; 

for (ang » 0; ang < ANGLE_STEPS/12 + AKGLE_STEPS/36 ; ang++) 
warped_angle [ang2 + ang] = angl + 1.25*ang; 
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angl = (3 *ANGriE_STBPS ) /4 + ANGIiE^STEPS/G ; 

ang2 = (3 *ANGIiE_STEPS ) /4 + ANGLE_STEPS/6 ' + ANGLE STEPS/36; 
for (ang « 0; ang < ANGLE__STEPS/l8 ; ang + + ) 
warped_angle [ang2 + ang] = aagl + l.5*ang; 

for (li = 0; li < I*_STEPS; /* initializes the UJT */ 

for (ai ^ 0; ai < A_STEPS; ai++) 
^ for (bi = 0; bi < B^STEPS; bi++) 

liUTCli] (ai) (bi] .c = (unsigned cha.r)0; 

LUTdi] (ai] Ibi] ,m = (unsigned cliar)0; 

LUTCliJ [ail [bi) .y = (unsigned char)0; 
inside_gainut [li] [aij [bi] = 0; 
concave [li] [ai] [bi] ^ 0; 



15 



25 



for (li r= lowL; li < L_STEPSr li++) /* sets the liUT CMYs */ 

for (ai = 0; ai < A_STEPS; ai++) /* inside printer gamut */ 

for (bi = 0; bi < B__STEPS; bi+ + ) 

20 t 

1 « i.nterval_l(li] ; 
a = interval_a [ai] ; 
b « interval_b [bi] ; 

if (ai !« CENTER_A || bi != CENTER__B) 

ang « angle (ai, bi) ; 
waorped^aing = warped_angle [ang] ; 
point. a ap- 
point, b = b; 

rotate (warped_ang - ang, point) ; 
a « rotated. a; 
b = rota ted. b; 

30 } 

fitd, a, b, &ci, &nii, &yi) ; 

if (ci >= 0 && ci <«= 64 mi >« 0 && mi <« 64 && yi >« 0 && yi 64 > 

convert ( ci , mi , y i ) ; 
IiUT[lil [ai] [bi] .c « cmy_vector . c ; 
35 LUT[li] [ai] [bi] .m * cray_vector . m ; 

IiUT[li] [ai] [bi] ,y «= cmy_vector .y; 
inside_garnut [li] [ai] [bi] = 1; 



ai CENTBR_A; /* resets the LUT CMYs */ 

bi = CENTER_B; /* at the L-axis */ 

for (li « 0; li < L^STEPS ; li++) 

ki = black [li] ; 
if (li >== lowL) 
45 inside_gainut [li] [ai] [bi] « i; 

IiUT(li] [ai] [bi] .c = LUT(li) [ai] [bi] .m « LUT(lil Cai] [bi] ,y 
- (unsigned char)lci; 

50 startii = L^STEPS; /* finds the printer gamut borders 
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25 



endh = -1; 

for (li = 0; li < L__STEPS; li++) 

startAdi] « A_STEPS; 

5 endACli] = -1; 

} 

for (li = 0; li < L__STEPS; li + + ) 

for (ai = 0; ai < A STEPS; ai-»-+) 
{ 

10 minB[li] [ai] = B_STEPS; 

maxB(li] [ai] = -l; 

for (li = 0; li < L STEPS; li++) 
{ 

15 for (ai = 0; ai < A STEPS; ai++) 

{ 

for (bi =0; bi < B STEPS; bi++) 
{ 

if (inside_gamut [li] [ai] [bi] === l) 

20 startL = li; 

break ; 

)' 

if (startL li) 
break; 

} 

if (startL == li) 
break; 

} 

for (li = L_STEPS - 1; li >= 0; li--) 
for (ai = 0; ai < A STEPS; ai-*-+) 

{ 

for (bi =0; bi < B STEPS; bi++) 

{ 

if (inside_gamut [li] [ai] [bi] == 1) 

endL = li; 
break ; 

if (endL li) 
break ; 

} 

if (endL li) 
break ; 

} 

45 

for (li ^ StartL; li <= endL; li++) 

for (ai = 0; ai < A_STEPS; ai+-H) 

for (bi =0; bi < B_STEPS ; bi++) 
50 { 

if (inside_gaxnut [li] [ai] [bi] =^ 1) 
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startAdi] s ai; 
break ; 



10 



15 



if (startACli] ai) 
break; 

) 

for (ai = A_STEPS - 1; ai >= 0; ai--) 

for (bi = 0; bi < B^STEPS; bi+-«-) 



{ 



if (inside_gamut [li] [ai] [bi] l) 



{ 



endAdi] = ai; 
break ; 



35 



if (endA[li] == ai) 
break ; 



20 



25 



30 



for (ai = startAUi]; ai <= endA[Xi] ; ai+-^) 
for (bi « 0; bi < B_STEPS ; bi+ + } 

if (inside_oramut [li] [ai] [bi] i) 



{ 



} 



raiuB[li] [ai] 
break; 



bi; 



} 

for (bi « B_STEPS - 1; bi >=» 0; bi--) 

if (inside_gainut: (li] [ai] [bi] i) 

rnaxB[li] [ai) «= bi; 
break ; 

, ) 



40 



for (li « startL; li endL; li++) 

for (ai = 0; ai < A STEPS; ai++) 

new_minB[li] [ai] = minB[liJ (aij ; 
^ new_iaaxB[li] [ai] = maxB[liI [ai] ; 



/* maJces Ctie printer gamut 
/* connected */ 



45 



50 



for (bi = minB[li] [CENTER_a] + l; bi < CENTER B; bi++) 



{ 



if {inside gamut [li] (CENTER_A) [bi] 
new_niinB[li] (CENTER A] ^ bi + 1; 



0) 



for (bi 
{ 



maxBdi] [CENTER_A] - 1 ; bi > CENTER_B; bi--) 
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if (insicie_gamut:[li] [CENTER_A] [bi] 0) 
new_maxB(li] [CENTER A] = bi - l; 



for 
{ 



(ai = CENTER_A + 1; ai <= endA[li] ; ai + + ) 



1] 
1]) 



for (bi = minBtli] [ai] + i; hi < inaxB[li] (ai] ; bi++) 
if (inside^gamut [li] [ai] [bi] == 0) 
new_minB[li] [ai] = bi + l; 
for (bi = maxB[li] [ai] - l; bi > minB[li] [ai] ; bi--) 
if (inside_gainut [li] [ai] [bi] 0) 
new_maxB[li] [ai] = bi - l; 
if (new_minB [li] [ai] > new_maxB [li] [ai] ) 

if (new_triinB[li] [ai] > new_maxB[li] [ai - 1] 
II maxB[li] [ai] < new__minB [li] [ai - 1]) 
new_minB[li] [ai] = minB[li] [ai] ; 
else if {minB[li) [ai] > new_maxB[li] [ai 
I I new__maxB [li] [ai] < new_ininB [li] [ai 
new_maxB[lil [ai] « maxB[li] [ai] ; 
else if (new_maxB [li] [ai] - minB[li] [ai] 
> maxB[li] [ai] - new_minB[li] [ai] ) 
new_rainB [li] [ai] = minB[li] [ai] ; 
else if (new_inaxB [li] [ai] minB[li] [ai] 
< maxBtli] [ai] - new_minB[li] [ai] ) 
new_niaxB(li] [ai] = inaxB[li] [ai] ; 
else if (minB[li] [ai] CENTER_B) 
new__minB [li] [ai] = minB[li] [ai] ; 
else 

new_inaxB[li] [ai] = maxB[li] [ai] ; 
if ( (iiew__minB [li] [ai] > new_nia:xB [li] [ai - 1] +1) 
II {new_inaxB[li] [ai] < new_minB [li] [ai - 1] - 1) ) 



{ 



for (aj ^ ai; aj <= endA[li] ; aj++) 

new_minB [ 1 i ] [a j ] « B_STEPS ; 
new_inaxB [li] [aj] = -1; 

break ; 
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for (ai = CENTER A - 1; ai >= startA[li] ; ai--) 
{ 

for (bi = minB[li] [ai] + 1; bi < maxB[li] [ail; bi++) 
if (inside_gamut: [li] [ai] [bi] " 0) 
new__minB [ 1 i ] [ai] = bi + 1; 
for (bi « maxB[li][ai] - 1; bi > itiinB [lil [ai] ; bi--) 
if (inside_gainut [li] [ai] [bi] 0) - 
new^maxB [li] [ai] = bi - l; 
if (new_minB(li] [ai] > new_inaxB [li] [ai] ) 

if (new^minB [li] [ai] > newjmaxB [li] [ai + 1] 
II majcB[li] [ai] < nGw_minB [li] [ai + 1]) 
new_minB [li] [ai] = minB(li] [ai] ; 
else if (minB(li] [ai] > new^maxB [li] [ai + 1] 
II new_maxB [li] [ai] < new__ininB [li] (ai + l] ) 
new_maxB [ 1 i 1 [ai] = maxB[li] [ai] ; 
else if (new_maLxB [li] [ai] - minB[li] [ai] 
> maxB[li] [ai] - new minB[li] [ai] ) 
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new_minB[li3 [ai] = minB[li] [ai) ; 
else if (new_maxBrii] [ai] - minB(lil tai] 
< maxBfli] (aij - new_minB [lij (ail) 
new maxBCli] (ai] = majcBLli] (aiJ ; 
- else if (minB[li] [ai] >« CENTER B) 
^0 new_minB[li] [aiJ = minB[li] [ai] ; 

else 

new_maxBfli] [ai] = majcBdi] [aiJ ; 

if ( (iiew_minB[li] [ai] > new_inaxB [Xi] (ai + i] + i) 
^ I I (nevr maxBdiJ [aij < new_minB [li J (ai -i- 1] - i) ) 

fo2r (aj = scartA[li] ; aj <= ai; aj++) 

new_minB[li] la j 3 = B^STEPS; 

^ new_inaxB(li] (aj ) = -in- 
break; 

20 J } 



25 



30 



for (ai « startAdi] ; ai <= endAdi] ; ai++) 

for <bi = minBdi] [ai] ; bi < new minB [li] [ai] ; bi++) 

inside_gainut [li] [ai] [bi] ^ 0; 
for (bi = new_ma3cB[li] [ai] + i; bi <= ma^cBdi] (ail; bi+ + ) 

xnside^gamut [li] [ai] [bi] « O; 
minBdi] [ai] = new_minB [li] [ai] ; 
^ maxBdi] [ai] = new_maxBdi] [ai] ; 

for (ai = CENTER_A + i; ai <= endAdi]; ai++) 

if (new_minB[liJ [ai] B_STEPS) 

for (aj = ai; aj endACli] ; aj++) 

for (bj = minB di] [aj] ; bj maxB[li] [aj] ; bj++) 

inside_gaimit [li] [aj] [bj ] O; 
endA[li] ^ b,± - i; 
brealc; 

>' 

for (ax = CENTER_A - 1; ai >= startAdi]; ai--) 
if (new_minB [li] [ai] B_STEPS) 

for (aj ^ ai; aj >= startAdi]; a j - - ) 

for (bj = minBdi] [aj]; bj maxBdi)[aj]; bj++) 
inside^gamut di] [aj] [bj] = 0; 
45 StartAdi] = ai + i; 

brealc; 

. > ' 



35 
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for (li = starti/; li <= endL; 11++) /* ma3ces the printer gamut 

^ /* radially convex */ 
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1 = interval__l [lil ; 

for (ang = 0; ang < ANGrjE_STEPS; ang+ + ) 

end^ai [angl = -1; 
end__bi[ang] = -i; 
warped_ang = warped_angle [ang] ; 
if (warped_ang < (ANGLE__STEPS) /8 || 
warped_ang >= {7*ANGLE STEPS) /8) 

{ 

ratio = tanpi ( (2 . 0*warped_ang) /ANGLE_STEPS ) ; 
for (ai = CENTER A + 1; ai < A STEPS; ai-«-+) 

{ 

a = interval_a [ai] ; 
b = ratio*a; 

bi = (int)((b - B_MIN) /DELTA_B) ; 
if (inside__gamut [li] [ai] [bi] o) 

start__ai [ang] ai; 
breaJc ; 

for (ai » start^ai [ang] + 1; ai < A_STEPS; ax++) 

a « interval_a [ai] ; 
b - ratio^a; 

bi - (int)((b - B_MIN) /DELTA_B) ; 
if (inside_gamut [li] [ai] [bi] == 1) 

end_ai [ang] = ai - 1; 
brea)c; 

else if (warped_ang >« (3*ANGLE_STEPS) /8 
&& warped_ang < (5*ANGLE_STEPS) /8) 

ratio « tanpi ( (2 . 0*warped__ang) / ANGLE_STEPS ) ; 

for (ai = CENTER A - 1; ai >= 0; ai--) 

{ 

a = interval_a [ai] ; 
b = ratio*a; 

bi - (int)((b - B_MIN) /DELTA_B) ; 
if (inside_gamut [li] [ai] [bi] == 0) 

start_ai [ang] = ai; 
brea)c; 

for (ai - start ai[ang] - 1; ai >= 0; ai--) 

{ 

a = inteirval_a [ai] ; 
b = ratio*a; 

bi - (int)((b - B_MIN) /DELTA_B) ; 
if (inside gamut [li] [ai] [bi] == 1) 

{ 

end_ai [ang] = ai + i; 
breaJc; 

) 

) 

) 



31 



EP 0 592146 A2 



else if (warped__ang >= (ANGI-E_STEPS) /8 
ScSc warped_ang < (3*ANGLE_STEPS) /8) 



{ 



if (warped_ang ANGLE_STEPS/4) 

ratio = 0; 
else 

ratio = l-0/tanpi( (2.0*warped_ang)/ANGLE_STEPS) 
for (bi = CENTER_B +1; bi < B_STEPS; bi++) 

b = interval_b [bi] ; 
a = ratio*b; 

ai - (int)<(a - A_MIN> /DEIiTA_A) ; 
if (inside_gainut (lil [ai] [bi] 0} 

start_bi [ang] = bi; 
breaJc ; 

for (bi = start_bi [smg] i; bi < B_STEPS; bi+ + ) 

b = iiiterval_b [bi] ; 
a « ratio*b; 

ai = (int)((a - A_MIN) /DELTA_A) ; 
if (inside_gamut [li] [ai] [bi] 1) 

end_bi[ang] = bi - 1; 
break ; 



} 

else 
{ 



if (warped_ang (3*ANGLE_STEPS) /4 ) 

ratio = 0 ; 
else 
{ 

^ ratio ^ 1 , 0/tanpi ( (2 . 0*warped_cLng) /ANGIjE_STEPS) 

for (bi = CENTER_B - 1; bi >= O; bi--) 

b = inte3rval_b [bi] ; 
a = ratio*b; 

ai - (int){{a - A_MIN) /DELTA_A) ; 
if (inside_gaiaut [li] [ai] [bi] ==0) 

start_bi [ang] = bi ; 
break ; 

> ' 

for (bi = start_bi [ang) - l; bi >= 0; bi--) 

b = interval_b [bi] ; 
a = ratio*b; 

ai = (int)((a - A_MIN) /DELTA_A) ; 
if (inside gamut [li] [ai] [bi] i) 

end_bi [ang] = bi + in- 
break; 
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} 

} 

5 for (ang = O; ang < ANGLE_STEPS; ang++) 

if (end_ai[ang] != -i || end_bi[ang] -1) 

warped__ang « warped_angle [ang] ; 
if {warped_ang < (ANGLE_STEPS) /8 
10 ^11 warped_ang (7*ANGLE_STEPS) /8) 

ratio = tanpi ( (2 • 0*warped_ang) /ANGLE_STEPS) ; 

for (ai = start_ai [ang] ; ai <= end_ai [ang] ; ai-i- + > 

a = interval__a[ai] ; 
15 b = rat:io*a; 

bi = (int)({b - B_MIN) /DELTA_B) ; 

down = intervcLl_b[bi] - 0,5*DELTA_B; 

up = interval_b [bi] + 0, 5*DELTA__B; 

left « ratio* (a - 0.5*DELTA_A) ; 

if (warped__ang < (ANGLE_STEPS) /8) 
20 { 

if (dovm - left > -0.001) 

first = bi - 1; , 
else 

first = bi; 

} ■ 

25 else 
{ 

if (left - up > -0,001) 

last = bi + 1; 
else 

last = bi; 

} 

right ^ ratio* (a + 0 . 5 *DELTA_A) ; 
i f ( warped_ang < ( ANGtiE_STEPS ) / 8 ) 

if (right - up > -0.001) 

last = bi + 1; 
else 

last = bi; 

) 

else 

{ 

if (dovm - right > -0,001) 

first = bi - 1; 
else 

first = bi; 

) 

for (bj = first; bj <= last; bj++) 

b = interval_b [bj ] ; 
if (inside_gamut [li] [ai] [bj] == 0) 

fitd, a, b, &ci, &mi, &yi) ; 
ci = suc)cin(ci) ; 
mi = suckindni); 
yi = suckin(yi) ; 
convert (ci, mi, yi) ; 
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LUTtli) Cai] [bj] .c = cmy_veccor.c; 
LUT[li] [ai] [b j 1 .m = cmy_vector .m; 
LOT[li) (ai] (bj] .y = ciny__vector .y; 
inside^gamut [li] (ai] (bj] « i; 
^ concave Cli] (ai] [bj ] = X; 

else if (warped__ang >= (3*ANGLE__STEPS) /8 
warped_ang < (5*ANGLiE__STEPS) /8) 



{ 



ratio «= tanpi ( (2 . 0*warped_ang) /ANGLE_STEPS) 
for (ai = starc_ai (ang] ; ai end_ai (ang] ; 

a =s iiitei-va.l_a [ail ; 
b = rat:io*a; 

bi = (int) ( {b - B_MIN) /DEIiTA_B) ; 
down - inte3rval_b[bi] - 0-5*DELTA_B; 
up = interval_b [bi] + 0 . 5*DEIjTA_B; 
left = ratio* {a - 0.5*DKIjTA__A) ; 
if (warped_ang {ANGIjE_STEPS) /2) 

if (dovm - left > -0,001) 

first = bi - 1; 
else 

first = bi; 

else 
{ 

if (left - up > -0.001) 

last = bi + 1; 
else 

last = bi; 

} 

right = ratio* (a + 0 .5*DEIjTA__A) ; 
if (warped_ang >« {ANGLE^STEPS) /2) 

if (right - up > -0,001) 

last = bi + 1; 
else 

last = bi; 

else 

{ 

if (down - right > -0.001) 

first bi - 1; 
else 

first = bi; 

} 

for (bj = first; bj <=- last; bj+ + ) 
b = interval_b (b j ] ; 

if (inside_gamut [li] [ai] (b j ] == 0) 

fit(l, a, b, &ci, fimi, &yi) ; 
ci = suckin(ci); 
mi = suckin(mi) ; 
yi = suckin(yi} ; 
convert (ci, mi, yi) ; 
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LUTCli] [ai] [bj] .c = cmy_vector .c.- 
LUTdi] (ai] [b j ] .m = cmy_vector .m; 
IiUT[li] [ai] [bjl.y = cmy_vector .y; 
inside_gaInut [li J [ai] (bj] - 1; 
concave[li] [ai] Cbj] = 1; 



else if (warped_ang >= (ANGLE_STEPS) /8 
&& warped__ang < (3*ANGLE_STEPS) /8) 

if (warped^ang ANGLE_STEPS/4) 

ratio = 0; 
else 

ratio = 1.0/tanpi( (2.0*warped_ang)/ANGLE_STEPS)* 
for (bi = start_bi [ang] ; bi <= end_bi[ang]; bi++) 

b = inteirval_b(bi] ; 
a = ratio*b; 

ai = (int)((a - A_MIN) /DELTA_A) ; 
down = interval^a [ai] - 0 . 5*DEIjTA_A; 
up = interval_a(ai] + 0 . 5 *DELTA__A; 
left = ratio* (b - 0 .5*DELTA_B) ; 
if <warped_ang <= (ANGLE_STEPS) /4) 

if (down - left > -0.001) 

first = ai - 1; 
else 
^ first = ai; 

else 
{ 

if (left - up > -0.001) 

last = ai + 1; 
else 

last = ai; 

} 

right = ratio* (b + 0. 5*DELTA__B) ; 
if (warped„ang <= (ANGLE_STEPS) /4) 

if (right - up > -0,001) 

last = ai + 1; 
else 

last = ai; 

} 

else 

{ 

if (down - right > -0.001) 

first = ai - 1; 
else 

first = ai; 

} 

for (aj = first; aj <= last; aj++) 

a = interval_a [a j ] ; 

if (inside_gainut [li] [aj] [bi] == 0) 

fitd, a, b, &ci, &mi, Scyi) ; 
ci = suclcin(ci) ; 
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mi = suckindni); 
yi = suckin(yi) ; 
convert ( ci , mi , yi ) ; 
LUT[li) (aj] [bi] ,c = any vector. c; 
LUTdi] [aj] (bil .m = cmy vector .m; 
liUTCli] [ajl Cbil .y = ciny__vector . y ; 
inside_gamut [11] [aj ] [bi] - X; 
concave (lil (ajj [bi] « 1; 



} 

Ise 

if (warped_ang ( 3 * ANGLE_STEPS ) /4) 

ratio = 0; 
else 

{ 

^ ratio = 1 - 0/ tanpi { {2 . 0*warped__ang) /ANGL.E_STEPS ) 

for (bi = start_bi [ang] ; bi end_bi [ang] ; bi--) 

b = interval_bEbi) ; 
a = ratio*b; 

ai = (int)((a - A_MIN) /DEIjTA_A) ; 
down = inteirval_a [ai] - 0 . 5*DEIjTA_A; 
up interval_a [ai] + 0.5*DEIjTA_A; 
left = ratio* (b - O . 5*DEIiTA_B) ; 
i f ( warped_ang <= (3 * ANGIiE_STEPS ) / 4 ) 

if (down - left > -0.001) 

first = ai - 1; 
else 

first = ai; 

else 
{ 

if (left - up > -0.001) 

last = ai + 1; 
else 

last = ai; 

} 

right = ratio* (b + 0 . 5*DELTA_B) ; 

if (warped_ang <= {3*ANGIiE_STEPS) /4) 

if (right - up > -0.001) 

last ai + 1; 
else 

last = ai; 

} 

else 

{ 

if (down - right > -0.001) 

first = ai - 1; 
else 

first = ai; 

} 

for (aj = first; aj last; aj+ + ) 

a = intervals [aj ] ; 
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if (inside_gainut: [li] [aj) [bi] =- 0) 

fit(l, a, b, &ci, &£ni, &yi) ; 

ci = suckin(ci) ; 

mi = suckin(mi) ; 

yi = suckin(yi) ; 

converted, mi, yi) ; 

LtJT[li] [ajj fbi] .c = cniy__vector.c; 

LUT[li] [aj] tbi] .m = czny__vector ,ni; 

LUT[li] (ajj [bi] ,y = cmy_vector , y ; 

inside_gamut [li] [aj] [bi] = 1; 

concave [li] [aj] [bi] = i; 



} 



} 



for (ai = startACli] ; ai <^ endA[lij ; ai++) 
{ 

minB[lil [ai] - B__STEPS; 
iwsLXB[li] [ai] = -1; 
for (bi = 0; bi < B STEPS; bi++) 

{. 

if (inside gamut [li] [ai] [bi] == 1) 

{ 

minB[li] [ai] - bi; 
brealc ; 

for (bi = B STEPS - 1; bi >=■ 0; bi--) 
{ 

if (inside__gainut [li] (ai] [bi] == 1) 

rnaxB[li] [ai] « bi; 
break; 

} 

for (ai = startAdi] ; ai <= endA[li] ; ai++) 

for (bi « minB(li] [ai] + 1; bi <= maxB[li] [ai] - 1; bi++) 

( 

if (inside^gamut [li] [ai] [bi] == 0) 
40 { • 

a ^ interval_a [ai] ; 
b = interval__b [bi] ; 
fit(l, a, b, ficci, &mi, &yi) ; 
ci = suckin(ci) ; 
mi = suckin(mi) ; 
45 yi = suckin{yi) ; 

convert (ci, mi, yi) ; 
LUT(li] [ai] [bi] -c cmy_vector .c; 
LUT(li] [ai] [bi] .m = cmy_vector .m; 
ljUT[li] [ai] [bi] .y = cmy_vector .y; 
inside_gamut [li] [ai] [bi] = i; 



} 
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) 

for (li « startL; li endL; /* extends the printer border 

for (ai = CENTER_A - 2 ; ai <- CBNTER_A + 2; ai^+) 

if (ai < startAtli] | | ai > endA(li3) 

minBtXi] [ai] = CENTER_B - 2; 
^ maxBdi] [ai] = CENTER_B + 2; 

else 
{ 

if (CENTER_B - 2 < minB[li] [ai] ) 

minB[lil [ai] = CENTER_B - 2; 
if {CENTER_B + 2 > maxB [li] [ai] ) 

maxBCli] [ai] = CENTER_B + 2; 

} 

if (CENTER_A - 2 < startAfli] ) 

startMli] = CENTER_A - 2; 
if (CSNTER_A + 2 > endA[li]) 

endA[li] = CENTER A + 2; 

} 

for (li = StartL; li <^ endL; li++) 

for (ai - 0; ai < startA[li] ; ai++) 

minB [li] [ai] = B^STEPS; 
maxBCli] [ai] = -i; 

} 

for (ai = endA[li] + i; ai < A_STEPS; ai4.+ ) 

minBCli] [ai] ~ B_STEPS; 
inaxB[li] [ai] =r -i; 

} ^ 

for (li « StartL; li <« endL; li++) 

startA[li] = startA[liJ - PRINTER_PLUS ; 
^ endA[li] = endA[li] + PRINTER_PLUS ; 

for (li = StartL; li <= endL; li++) 

for (ai « startA[li]; ai <« endA[li] ; ai++) 

new_minB[li] [ai] = minB[li] [ai] - PRINTER_PLUS ; 
new__niaxBfli] [ai] = maxB[li] [ai] + PRINTER_PLUS ; 
for (aj « ai - PRINTER_PLUS ; aj <= ai + PRINTER_PLUS ; aj++) 

if (aj >= 0 && aj < A_STEPS) 

new_minB [li] [ai] 

= min2 {new_minB [li] [ai] , ininB[li] (aj ] - PRINTER PLUS) ; 
new_maxB[li] (ai] ~ 

^ = max2 (new_maxB[li] [ai] , maxB [li] (a j ] + PRINTER_PLUS ) ; ' ' 
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for (ai ^ st:artA[li]; ai endAfli] ; ai+^) 

minBdi] [ai] « new_minB [li] [ai] ; 
maxBtli] [ai] = new_maxB[li] [ai] ; 

} 

for (li = startL; li <= endL; li++) 

m_,startA[li] = m_startAUi] - MONITOR PLUS; 
^ m_endA[li] - m^endAtli] + MONITOR_PLUS; 

for (li ^ StartL; li <= endL; li4.+) 

for (ai = m^startAfli] ; ai <= in_endAtli]; ai^+) 

new__minB[li] [ai] « m_minB[li3 [ai] - monitor PLUS; 

new_,inaxB[li] [ai] « m^inaxB[li] [ai] + MONITOR~PLUS 

for (aj ai - MONITOR_PLUS ; aj ai + MONITOR_PLUS ; aj+^) 

if (aj >= 0 && aj < A_STEPS) 

new_minB[li] [ai] 

new_^[li7TSr^^^^ rn_minB[li] [aj] - MONITOR_PLUS ) 

^ « inax2(new_maxB[li] (ai] , m_maxB[li] [aj] + MONITOR_PLUS ) 

for (ai = in_startACli] ; ai <= m_endA(li]; ai++) 

in_minB[li] [ai] = new_minB(li] [ai] ; 
^ m_inaxB[li] [ai] = new__inaxB [li] [ai] i 

for (li = StartL; li <= endL; li++) 
for (ai = 0; ai < A_STEPS; ai++) 

^^[^•1 {^'1 " min2(minB[li] [ai] , m_niinB[li] [ai] ) ; 
^ majcB[li] [ai] = inax2 (maxB [li] [ai] , m~maxB [li] [ai] ) ; 

for (ai = 0; ai < A_STEPS; ai++) 

if (minB[li] [ai] < B^STEPS && maxB [li] [ai] > -i) 

startA[li] = ai; 
break; 

for (ai = A_STEPS - 1; ai >= 0; ai--) 

if (minBLli] [ai] < B^STEPS && maxB [li] [ai] > -i) 

endAdi] = ai; 
breaJc; 

for (li = StartL; li <= endL; li++) 

for (ang 0; ang < ANGLE_STEPS; ang++) 
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if {{ang < ANGLE_STEPS/8 || ang (7*ANGLE_STEPS) /8 ) || 

^ (ang >= (3*ANGIiE_STEPS) /8 ang < <5*ANGLE_STEPS) /8 ) ) 

ratio = tanpi( (2.0*ang)/ANGLE_STEPS) ; 

if (ang < ANGLE_STEPS/4 || ang >= (3*ANGLE_STEPS) /4) 

for (ai = A_STEPS - l; ai >= CENTER_A; ai--> 

a = interval_a (ai] ; 
b = ratio*a; 

bi = (int:)((b - B_MIN) /DELTA_B) ; 

if (bi >= minB[li) [ai] &£e bi <= n\axB[li] [ai] ) 

outer^border [li] [ang] ,a = a; 
outer_border [lij [ang] .b = b; 
break; 

else 

{ 

for (ai = 0; ai <« CENTER_A; ai + + ) 

a = interval_a [ai] ; 
b = ratio*a; 

bi = (int)((b - B_MIN) /DELTA_B) ; 

if (bi >= minBdi] [ai] && bi <^ maxB[li] [ai] ) 

outer_border [li] [ang] .a = a; 
outer_border [li] [ang] .b = b; 
break; 

y' 

else 
{ 

if (smg «= ANGIiE_STEPS/4 | | ang (3 *ANGLE_STEPS) /4) 

ratio = 0 ; 
else 

ratio = 1.0/tanpi( (2.0*ang)/ANGLE_STEPS> ; 
if (ang < ANGIiE_STEPS/2 ) 

for (bi = B_STEPS - 1; bi >= CENTER_B; hi--) 

b ~ interval_b (bi] ; 
a = ratio*b; 

ai = (int)((a - A_MIN) /DELTA_A) ; 

if (bi >= minBdi] [ai] && bi <= maxBdi] (ail) 

outer_border [li] [ang] .a = a; 
outer_border [li] [ang] .b = b; 
break; 

else 

{ 

for (bi = O; bi <= CENTER B; bi++) 
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{ 

5 b = interval_b [bi] ; 

a = ratio*b; 

ai ^ (int:)(<a - A_MIN) /DELTA_A) ; 

if (bi >= minBCli] [ai) && bi <= maxB[li) [ai] ) 

outer_border [li] [angj .a = a; 
outer__border [li] [angJ .b = b; 
10 break; 



15 

for fang = 0; ang < ANGLE_STEPS; ang++) /* makes vertical outer */ 

{ /* border profiles monotonic 

if (aJig < ANGLE_STEPS/4 || ang > <3*ANGliE_STEPS) /4 ) 
aa = 1; 

if (ang > ANGLE_STEPS/4 && ang < (3*ANGLE_STEPS) /4 ) 
aa - -1; 

if (ang ANGLE_STEPS/4 || ang (3*ANGLE_STEPS) /4) 
aa « 0; 

if (ang > 0 && ang < ANGLE_STEPS/2) 
bb = 1; 

if (ang > ANGLE_STEPS/2 && ang < ANGLE^STEPS) 
25 bb = -1; 

if (ang 0 1 | ang ANGIiE„STEPS/2) 

bb - 0; 
11 = saturation_l [ang] ; 
a = outer_border [11] [ang] .a; 
b = outer_border [11] [ang] .b; 
for (li = 11 + 1; li <= endli; li++) 

{ 

if (aa*outer_border [11] [ang] .a > aa*a 
I I bb*outer_border[ll] (ang) .b > bb*b) 

outer_border [11] [ang] .a = a; 
outer_border [11] [ang] .b = b; 

35 } 

a - outer_border [11] [ang] .a; 
b = outer__border [11] [ang] .b; 

a ^ oucer_border (111 (ang] .a; 
b - outer_border [11] (ang] .b; 
for (li « 11 - 1; li >= startL; li--) 
{ 

if (aa*outer_border [11] (ang] .a > aa*a 
|j bb*outerborder [11] [ang] .b > bb*b) 

outer_border [11] [ang] .a = a; 
outer_border [11] (ang) ,b = b; 

45 } . 

a = outer_border (11] (ang] .a; 
b = cuter border [11] [ang] .b; 

) 
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20 



for (ang » 0; ang < ANGLE_STEPS; angi. + ) /* sets the saturation */ 

5 ^ saturation = -X; border LABs and CMYs */ 

warped_ang « wajrped_angle Cetng] ; 

if <warped_ang < AKGLE_STEPS/a \ \ warped_ang >= (7*angle_STEPS) /8 ) 

ratio = tanpK (2.0*warped_ang) /ANGLE STEPS); 
point. a = interval__a fA_STEPS - 1] ; 
10 ^ point, b = ratio*point ,a; 

else if (warped^ang >« ( 3 * ANGLE_STEPS ) / 8 && warped_ang < (5*ANGLE_STEPS) 

ratio = tanpi ( (2 .0*warpeci_ang) /ANGLE STEPS)- 
point. a » interval_a[0 3 ; ~ 
^ point.b = ratio*point ,a; 

else if (warped_ang >« ANGLE_STEPS/8 && warped_ang < (3*ANGLE_STEPS) /8) 

if (warped ang «= ANGLE_STEPS/4 ) 

ratio = 0; 
else 

ratio = l-0/tanpi( (2-0*warped__ang)/ANGLE STEPS); 
point.b = interval__b (B_STEPS -13; ~ 
^ point. a = ratio*point.b; 

else 
{ 

25 (warped_ang (3*ANGIiE_STEPS) /4) 

ratio « 0; 
else 

ratio l.O/tanpi( (2,0*warped ang) /ANGLE STEPS); 
point. b = interval_b[0] ; " 
^ point. a = ratio*point .b; 

3^ for (li - startL; li encLL; li++) 

Im = intervai_l (li] ; 

intersection (Im, ang, point, 0. &1, &a, &b) ; 
new_saturation « a*a + b*b; 
if (new_saturation > saturation) 
35 { 

11 - li; 

^ saturation new_saturation; 

saturation = -i; 
^ if (11 == StartL) 

Im - interval_l cm ; 
else 

Im =» interval_l [11 - i] ; 
delta_l = DELTA_L/64; 
if (11 StartL II 11 == endL) 

last = 64; 
^5 else 

last = 128; 
for (index = 1; index < 12 8; index+ + ) 

intersection dm + index*delta_l , ang, point, 0, &1, &a, &b) ; 
new_saturation = a*a + b*b; 
50 if (new_saturation > saturation) 
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10 



15 



) 



} 



In = 1; 

an = a; 
bn » b; 

saturation = new saturation; 



saturation_l [ang] = In; 
saturation_border (angj , a = an ; 
saturation_border [ang] .b = bn; 
fit (In, an, bn, &ci, &nii, &yi) ; 
final_ci = suckin(ci) 
final_mi = suckin(mi) 
finally i « suckin{yi) 
convert (final_ci, final_mi, final_yi) ; 
saturation_border_color [ang] . c cmy_vector-c; 
saturation_border_color (ang] .m = cmyvector .m; 
saturation_border_color lang] .y « cmy_vector.y; 



20 
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for (li » startL; li <« endL; li++) 
{ 

Im ss interval_l (lij ; 

for (ang == 0; ang < ANGLE_STBPS; ang++) 



/* sets the border CMYs */ 



{ 



point* a outer__border [li] tang] .a; 
point -b = outer_border tli] (ang] .b; 
warped__ang « warped_angle (ajig] ; 
rotate (warped_ang - ang, point); 
point. a « rotated, a; 
point. b = rotated. b; 
radm = radius (point) ; 

rad = radius (saturacion_border [ang] ) ; 
In « saturation_l [ang] ; 
if (in < Im) 

sign 1; 
else if (In > Im) 

sign = -1; 
else 

sign = 0; 

if (radm > rad && (sign*(lm - In)) /(radm - rad) < SPACE_RATIO) 

bordercolor [li] (ang] .c = saturation_border_color [ang] ,c; 
border_color [li] (ang] .m ^ saturation_border_color [ang] .m; 
border_color(li] (ajig] .y « saturation__border__CQlor (auag] .y ; 

else 

{ 

intersection dm, ang, point, sign, sd, Sa, &b) ; 

if (a 0 && b -= 0) 

{ 

for (Ij = StartL; Ij <= endL; lj++) 

if (interval_l [Ij] + (interval_l [1 j + 1] 
breaJc; 

} 

if (Ij == endL + i) 

)cj = 0; 
else 



interval_l CI j 1 )/2.0 
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kj » biackClj] ; 

bord^^-^''^''''f^'^^ ^^9^ - (unsigned char) kj 
H^^H -''''^•''''^^^ -"^ « (unsigned char) )cj 

^ border_color[li} [ang] .y - (unsigned char) kj , 

else 
{ 

fitCl, a, b, &:ci, &mi, iyi) ; 
fxnal ci t= suckin<ci) ; 
final_mi = suckin(mi) ; 
f inal_yi « suckin(yi) ; 

convert (final_ci, final_nii, f inalyi) ; 
border_colorlli] [ang] .c = cmy__vector.c; 
border_colorrii) (ang) .m « cmy_vector .m; 
border_color[li] (ang] .y - cmy^vector.y ; 



for (li = startL; li endL; li+^) /* converges the LOT into 

^ im ^ interval,! [1x3; ^* ^""^ ''^^^^^ 

for (ai « startAdi); ai <= endAtli] ? ai+ + ) 

for (bi = minBdi] [ai] ; bi maxB(li] [ai] ; bi+ + ) 

point. a « interval_a[ai] ; 

point, b = interval_b[bi] ; 

saturation = radius (point ) ; 

if (ai CE^FTER__A 11 bi 1= CENTER_B) 

ang angle (ai, bi) ; 

warped_ang = warped_eLngle [ang] ; 

rotate {warpedang - ang^ point); 

point. a = rotated, a; 

point. b « rota ted. b; 

if {inside_gamut [li] [ai] Cbi] o) 

radm = radius (point ) ; 

rad « radius (saturation_border [ang] ) ; 
In = saturation_l [ang] ; 
if (In < Im) 

sign « 1; 
else if (In > Im) 

sign = -1; 
else 

sign = 0; 

if (radm > rad && (sign*(lm - ln))/(radm - rad) < SPACE_RATIO) 

liUTUi] (ai] [bi] .c - saturation_border_color [ang] .c; 
i'J!! . ^ l^^l ^^^^ saturation_border_color[ang] .m; 

liUTili] [ai] [bi] .y « saturation border color(anQl v 
^ transition tli) [ai] [bi] = i; ~ " ^ -J^ 

else 
{ 

intersection dm, ang, point, sign, il . &a, &b) ; 
if (a 0 && b 0) 
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10 



15 



20 



{Ij = startL; Ij 



for 
{ 

if (interval_l (Ijl 
break; 

) 

if (Ij endL + l) 

kj = 0; 
else 

kj = black [Ij] ; 
LUT[li] [ai] [bi] .c = 
LUT[lil [ail [bi] ,m = 
LTJT(li) [ai] [bi] .y - 
transition tli] [ai] [bi] 

} 

else 

{ 



<= endL; lj++) 
+ DELTA L/2.0 > 1) 



(unsigned char) kj ; 
( uns igned char ) kj ; 
(unsigned char) kj ; 
= 1; 



fitd, a, b, &ci, &nii, &yi) ; 
final_ci = sucfcin(ci>; 
f inal__mi. = suckin(mi); 
final_yi = suckin(yi) ; 

convert (f inal_ci, f inal_iai , f inal_yi) ; 
LUT[li3 [ai] [bi] .c « cmy_vector . c? 
LUTdi] [ail tbi] .m = ciny_vector .m; 
LUT[li] [ai] [bi] .y « cmy_vector-y ; 
transition [li] [ai] [bi] « 1; 
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for 
{ 



(index = 0; index < 2; index++) 



close_a [index] 
close b [index] 



CENTER_A + 1; 
CENTER B + index; 



close_a [2 
close_b[2 
close_a [4 
close_b [4 
close_a [6 
close_b[6 

} 

for {index = 

{ 



index] ~ CENTER_A - index; 

index) = CENTER_B + 1; 

index] ^ CENTER_A - 1; 

index] = CENTER_B - index; 

index] = CENTER_A + index; 

index] « CENTER_B - 1; 

0; index < 4; index+-*-) 



far_a [index] = CENTER_A + 2; 
far_bC index] = CENTER_B - 1 + index; 
far__a[4 + index) = CENTER_A + l - index; 
far_b[4 + index] = CENTER_B + 2; 
far_a[8 + index] = CENTER_A - 2? 
far_b[8 + index] = CENTER_B + 1 - index; 
far_a[l2 + index] « CENTER_A - 1 ■♦• index; 
far b[12 + index] = CENTER_B - 2; 



/* makes- the colors near */ 
/* the gray axis more gray 



for (index = 0; index < 5; index++) 



{ 



50 



border_a [index] = CENTER_A + 3; 
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border 
border 
border 
border 
border 
border 
border 

} 

border_a 
border__b 
border_a 
border_b 
border__a 
border__b 
border_a 
border b 



b (index) ^ CENTER_B - 2 + index; 
_a[6 + index) « CENTBR_A + 2 - index; 
'__bC6 -t- index) = CENTER_B + 3; 
■_a[12 + index] « CBNTER_A - 3; 
•_b(X2 + index) = CENTER_B + 2 - 
• aCie + index] = CENTER_a - 2 + 
b[18 + index) = CENTER B - 3; 



index; 
index; 



C5] = CENTER_A + 3; 
C5) = CENTER_B + 3; 
(113 = CENTER^A - 3 
(11] = CENTER_B + 3 
(17] « CENTER_A - 3 
(17] = CENTER_B - 3 
(23] = CENTER_A + 3 
(23] = CENTER B - 3 



20 



for 
{ 



(li 



startL; li endL; li++) 



1 = interval__l (li] ; 
for (index = 0; index < 8; index++) 
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ai = border_a [3*index +2]; 
bi « bor<ier_b t3*index + 2] ; 
aa = far__a[2*index + U ; 
bb = far_b(2*index + 1]; 
ang = angle (aa, bb) ; 
if (inside^gamut: [li] [ai] [bi] 



{ 

hue [2* index 
hue (2* index 
hue [2* index 

} 

else 

{ 

hue ( 2 * index 
hue [ 2 * index 
hue (2*index 



1 I I transition (li) [ail [bi] l) 



+ 1] .c = LOTdi] [ai] [bi] .c; 
+ 1] .m = LUTdiJ [ai] [bi] .m; 
+ 1] .y « LUTdiJ [ail Cbi] .y; 



1] .c = border_color (li) (ang] 
1] -m = border__color [li] (ang] 
1] .y « border_color [li] (ang] 



c; 
m; 

y; 



} 

k]c ^ min3( <int)hue(2*index + l].c, hue(2*index + 1] -m, 

hue(2*index + l] .y) ; 
if (ang > n\agenta_ang | | ang < yellow_ang) 

hue(2*index + ll.c - (unsigned chax)k)c; 
if (ang > yellow_ahg && ang < cyan^ang) 

hue(2*index + 1] ,m = (unsigned char))c)c; 
if (ang > cyan_ang && ang < magenta_ang) 
hue(2*index + l) .y = (unsigned char)lck; 

index < 8; index++) 



} 

for (index = 0; 
{ 

ai = border_a (3*index] ; 
bi = border_b{3*index] ; 
aj = bordera (3*index + 
bj = border_b (3*index + 
aa = f ar_a [2*index] 
bb - far_b [2*index] 
ang - angle (aa, bb) 



1] ; 
1] ; 
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average2(li, ai , bi , aj , b j , 
hue (2* index] . c « averaged. c; 



0) 
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hue[2*index] .m « averaged. m; 
hue (2* index] .y = averaged- y; 

kk min3 ( (int)hue(2*index] .c, hue [2* index] ,m, hue (2* index] .y) ; 
if (ang > magenta_ang | | ang < yellow_ang) 

hueC2*index] .c « (unsigned ch^r)kk:; 
if (ang > yellow_ang && ang < cyan_ang) 

hue [2* index] .m = (unsigned char) kit; 
if (ang > cyan_ang ang < magenca_ang) 

hue [2* index) .y ^ (unsigned char)kk; 

10 } 
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ki = min3 ( (inc) LUTCli] [CENTER^A] [CENTER_B] .c, 

(int)IiUT(li] CCEWTER_A] [CENTBR_B] .m, (int)LUT(li] [CENTER_A] [CENTER_B] .y 
for (index = 0; index < 8; index++) 

15 I 

cc « (int)huet2*index + 1] -c; 
mm = (int) hue [2-^ index + 1] .m; 
yy « (int)hue[2*index + 1] .y; 

Ick = minS ( (int)hue(2*index + l].c, (int) hue [2* index + 1] .m, 

(int) hue (2*index + l] .y) ; 
ci = (int) (ki + (1 . 0/3 . 0) *( (double) (cc - ki) ) ) 
mi = (int) (ki + (1 , 0/3 . 0) *( (double) (mm - ki))) 
cj - (int) (ki + (2,0/3.0) *( (double) (cc - ki))) 
mj = (int) (ki + (2 . 0/3 . 0) * ( (double) (mm - ki))) 
if (ki <= 1 && (ci == cj II mi ==mj) && (cc != 0 || mm !«= 0)) 

yi « (int) (ki + (2 . 0/3 , 0) *( (double) (yy - ki) ) + 0.5); 
else 

yi = (int) (ki + (1 . 0/3 . 0) *( (double) (yy - ki) ) + 0.5); 
if (ci == mi && mi == yi && ci != 0) 

if (cc == kk) 

ci - - ; 
if (mm == )ck) 
30 mi--; 

if (yy == kk) 
yi--; 

} 

kj « min3(ci, mi, yi) ; 
if (ki <= 1 && kj > kk) 

35 { 

if (cc == kk) 

ci = kk; 

if (mm === kk) 

mi - )ck; 

if (yy == kk) 

yi = kk; 



if (ci == 0 && mi ^= 0 && (ki != 0 | | index 2)) 

if (cc > mm) 

ci + + ; 
else if (cc < mm) 

mi + + ; 
else 

{ 

ci + + ; 
mi + + ; 

) 
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ai « close__a [index) ; 
bi = close_b [index] ; 

LUTUi] [ai] (bij .c ^ (unsigned char)ci; 

IiUTlli] [ail (bi] .m = (unsigned char) mi; 

LUT[li] [ai] [bi) .y = (unsigned char)yi; 
inside_gamut [li] [ai] tbi] = 1; 



20 



10 for (index « O; index < 16; index+ + ) 

cc = (int) hue [index] .c; 
mm = (int) hue [index] .m; 
yy = { int ) hue [ index) . y ; 

kk - min3 ( (int) hue [index) .c, (int) hue [index] .m, ( int ) hue [index] .v) - 
15 = (mt) (}ci + (2.0/3.0) *( (double) (cc - ki) ) ) ; ' 

mi = (int) (ki + (2 . 0/3 , 0) * ( (double) (mm - ki) ) ) ; 
yx « (int) (lei + (2 , 0/3 , 0) *( (double) (yy - ki) ) + 0.5); 
if (ci mi && mi =- yi ci != 0) 

if (cc kk) 

ci--; 
if (mm == Ick) 

mi - - ; 
if (yy ]ck) 

kj « min3(ci, mi, yi) ; 
25 if (ki <= 2 && kj > kk) 

{ 

if (cc )ck) 

ci = kk; 
if (mm == ]ck) 
mi = )ck; 

if (yy == )ck) 
^ yi . - kk ; 

if (ci 0 && mi ~ 0 && (ki S= 0 | | (index != 5 && index !« 6))) 
if (cc > mm) 

Ci + +; 

else if (cc < mm) 

mi++ ; 
else 

{ 

ci + +; 
mi + + ; 

40 } 

} 

ai => far_a [index] ; 
bi = far_b [index] ; 

LUT[li] [ai] [bi] .c = (unsigned char)ci; 
LUTCli] [ai] [bi] .m (unsigned char)mi; 
45 LUT[li] [ai] [bi] .y = (unsigned char)yi; 

inside_gamut [li] [ai] [bi] = i; 
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for (index = 0; index < 5; index++) /* recomputes colors in */ 
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10 



^ For- n ■! - c^f^^fT . /* concave areas */ 

for [11 ^ startL; li <^ endL; li-f. + ) ' 

^ ^^^wT startAEli]; ai <= endAdi] ; ai++) 
for (bx « miiiB[li] [ai]; bi cnaxB [li] [ai] ; bi++) 



[ 



if (concave [li] [ai] [bi] == i) 

average4(li, ai, bi, i) ; 
LUT(li) [ai] [bi] .0 = averaged. c; 
LUT[li] [ai] [bi] .m = averaged. m; 
LUT[li] [ai] [bi] .y « averaged. y; 



15 



20 



for (li ^ StartL; li <« endL; li+ + > 

starts [li] = minB [li] [CENTER A] ; 
endBdi] « inaxB[liJ [CENTER_A] ; 
for (ai = startA[li]; ai <= endA[li] ; ai++) 



{ 



if (minB[li] [ai] < startB[li]) 
startB[li] « minBCli] [ai] ; 

if (inaxB[li] (ai] >endB[li]) 
endB[li) = maxB [li] [ai] ; 



/* computes startB, endB */ 
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for (li 
{ 



StartL; li <= endL; li++) 



for (ai = startAdi]; ai <= endA[li] ; ai++) 



( 



/* fills con^lete LUT */ 
/*- rectcuigles */ 



startB[li3; bi < minB [li] [ai] ; bi++) 



for (bi 
{ 

ang = angle (ai, bi) ; 

HJT[li] [ai] [bi] .c := border_color [li] (ang] .c; 
LtJT[li] [ai] [bi] .m = border_color [li] [ang] .m; 
^ LUTCli] [ai] [bi] .y = border_color[li] [ang] .y; 

for (bi - inaxB(li](ai] + i; bi <= endBdi]; bi++) 
ang « angle (ai, bi) ; 

LUT[li] [ai] [bi] .c « border_color (li) [ang] .c; 
LUTdi] [ai] [bi] .m = border_color [li] [ang] .m; 
LUTdi] [ai] [bi] ,y ^ border_color [li] [ang] .y; 



45 



for (li = 0; li < StartL; li++) 



StartAdi] « endA[li] ;= CENTER A; 
starts [li] = endBdi] = CENTER_B; 



/* sets startA, endA, startB, endB * 
/* and border CMYs for low Ls */ 



for (11 = 0; li < StartL; li++) 

for (ang = 0; ang < ANGLE__STEPS ; ang++) 
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10 



15 



20 



\^ll±V(r^^l^'^.^^^^ £CENTER_A3 (CENTER B] .c, 
\ ^^^^-^^ [CENTER B]Tm, - ^ 

Hr.i;!?Ii [CBNTER-B] ,y); 



} 

fwrite(&startL, sizeoffint) i /* . 

lab_borders jptr) ■ ^ writes the lab__borders */ 

ror til « 0; li < L^STEPS; + -i* " 

^ rwriteCtendAdi] , Sizeof(int), i, lab_borders_ptr) ; 
for (li = 0; li < I._STEPS; li++) 

l^itSj^n^^'pf f^^^ • {i-it^) - 1. lab borders_ptr) - 

^ fwrite(&eadB(li] , sizeof(int), i, lab_b5rders_pt?) ; 

for (li = 0; li < L STEPS- li++) /. 

'"for'kT -"^""^f^i^.- ai ci :ndA[li]; ai.., /- """^"^ 

for (bl « startBdi]; bi <« t^nrfRnsi. -k-j ^ 
f write («,trr Cli] tai/ [bii , ''sizeol^^^f i , i^tjtr, ; 

25 for (li « 0; li < L STEPS- li-*.-*.^ /* 

for (ang = O; ang-< ANgLe S^EPS; ang..) ^ "''^''^^ ^^'^ border_colors 

fwrxte(a3order_colortli][Sag]. sizelf (cmy) , i. border_colors_ptr) ; 
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35 Claims 



1. A color printer look-up table for providing color primary values corresponding to colors in device-indepen- 
dent color space, the color printer table having smoothly varying color values for colors outside the color 
printer gamut. 

40 

2. A color printer look-up table according to Claim 1, wherein colors outside the color printer gamut have 
monotonically increasing lightness from the darkest value in the printer table to the lightest value in the 
printer table. 

^5 3. A color printer look-up table according to Claim 1 or 2, wherein the color primary values include color pri- 
mary values for colors inside the color printer gamut and for colors outside the color printer gamut. 

4. A color printer look-up table according to Claim 3, wherein the color primary values include color primary 
values for colors within a color monitor gamut 

50 

5- A color printer look-up table according to Qalm 3 or 4, wherein for colors inside the color printer gamut, 
the color primary values are provided In accurate correspondence to the associated colors. 

6. A color printer look-up table according to Claim 5, wherein the color primary values outside the color printer 
gamut are provided so as to preserve hue. increase saturation, and change lightness. 

7. A color printer driver comprising: 

a printer table for providing color primary values corresponding to colors within the color printer 
gamut and for colors outside the color printer gamut the color primary values varying smoothly for colors 
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outside the color printer gamut; 

means for accepting a command for printing a designated color; and 

control means for extracting color primary values corresponding to the designated color from the 
printer table. 

8. A color printer look-up table according to Claim 7, wherein colors outside the color printer gamut have 
monotonically increasing lightness from the darkest value in the printer table to the lightest value in the 
printer table. 

9. A color printer driver according to Claim 7 or 8, wherein the printer table includes color primary values for 
colors within a color monitor gamut. 

10. A color printer driver according to Claim 7, 8, 9 or 10 wherein for colors within the color printer gamut, the 
color primary values in the printer table are provided in accurate correspondence to the associated colors, 

11. A color printer driver according to any of Claims 7 to 10, wherein the color primary values for colore in the 
printer table outside the printer gamut are provided so as to preserve hue, increase saturation, and change 
lightness. 

12. A color printer driver according to any of Daims 7 to 11, further comprising a border table for providing 
color primary values for colors outside the printer table, wherein said control means selects the printer 
table or the border table based on the designated color and extracts color primary values from the selected 
one of the printer table and the bonder table. 

13. A color printer driver according to Claim 12, wherein the border table is arranged in a wheel-like arrange- 
ment of cells centered on the lightness axis. 

14. A color printer driver according to Claim 12 or 13, wherein the color primary values for colors in the bonder 
table outside the printer gamut are provided so as to preserve hue, increase saturation, and change light- 
ness. 

15. A method for building a printer table comprising: 

determining a color printer gamut edge; 

mapping color primary values into the printer table for colors within the color printer gamut edge; 

and 

mapping transition colors into the printer table between the printer gamut edge and the edge of 
the printer table by constant angle extension to the edge of the printer gamut edge so as to increase color 
saturation. 

16. A method according to Claim 15, further comprising the steps of mapping all points in a wedge subtended 
by the constant angle from the maximum saturation point of the printer gamut edge onto the maximum 
saturation point. 

17. A method according to Claim 15 or 16, wherein the constant angle is substantially about 15^ 

18. A method according to Claim 15, 16, 17 or 18, wherein the color primary values for transition colors in 
the printer table are provided for colors within a color monitor gamut. 

19. A method according to any of Claims 1 5 to 1 8, wherein for colors within the color printer gamut, the color 
primary values are provided in accurate correspondence to the associated colors. 

20. A method according to Claim 1 9, wherein the color primary values outside the printer gamut are provided 
so as to preserve hue, increase saturation, and change lightness. 

21. A method according to any of Claims 15 to 20, further comprising the step of building a border table. 

22. A method according to Claim 21 , wherein the color primary values for colors in the border table outside 
the printer gamut are provided so as to preserve hue, increase saturation, and change lightness. 

23. A method of color image repnaduction wherein not all desired colore are available via a desired output de- 
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vice, and wherein unavailable colors are mapped to respective availabi colors so as to preserve mono- 
tonic increases in lightness of the desir d color. 

A method of color image reproduction wherein not all desired colors are available via a desired output de- 
vice, and wherein unavailable colors are mapped to respectiv. available colors so as to preserve an angle 
in color space between the desired color and a plane of constant luminance. 

An image signal or recording generated using a method or apparatus according to any preceding claim. 
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@ Color reproduction method and apparatus. 

@ Method and apparatus for color printing ac- 
cording to a printer table having high color 
smoothness for out-of-gamut colors. A color 
printer gamut edge is first determined and color 
primary values for colors within the printer 
gamut are calculated and inserted into the 
printer table. For transition colors, namely 
those colors outside the printer gamut edge but 
still within the printer table, color primary 
values are calculated by selecting a color from 
the printer gamut edge at a point which lies at a 
constant angle from the transition color in 
question. The same constant angle is used for 
each and every one of the transition colors. A 
border table is provided for colors outside of 
the printer table ; the color primary values for 
the border table are selected in the same man- 
ner as for the transition colors. According to the 
invention, it is possible to provide smooth color 
transitions for out-of-gamut colors, and in par- 
ticular it is possible to provide monotonically 
increasing lightness in both the printer table 
and the border table and thereby avoid undesir- 
able bands of darker colors within lighter col- 
ors. 
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